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Abstract 

We present the first unquenched lattice-QCD calculation of the hadronic form factors for the 
exclusive decay B —)• Dil^ at nonzero recoil. We carry out numerical simulations on fourteen 
ensembles of gauge-field configurations generated with 2-\-l flavors of asqtad-improved staggered 
sea quarks. The ensembles encompass a wide range of lattice spacings (approximately 0.045 to 
0.12 fm) and ratios of light (up and down) to strange sea-quark masses ranging from 0.05 to 0.4. 
For the b and c valence quarks we use improved Wilson fermions with the Fermilab interpretation, 
while for the light valence quarks we use asqtad-improved staggered fermions. We extrapolate our 
results to the physical point using rooted staggered heavy-light meson chiral perturbation theory. 
We then parameterize the form factors and extend them to the full kinematic range using model- 
independent functions based on analyticity and unitarity. We present our final results for f+{q‘^) 
and fo{q‘^), including statistical and systematic errors, as coefficients of a series in the variable z and 
the covariance matrix between these coefficients. We then fit the lattice form-factor data jointly 
with the experimentally measured differential decay rate from BaBar to determine the CKM matrix 
element, \Vcb\ = (39.6 ± l.TgcD+exp i 0.2 qed) x 10“^. As a byproduct of the joint fit we obtain the 
form factors with improved precision at large recoil. Finally, we use them to update our calculation 
of the ratio R{D) in the Standard Model, which yields R{D) = 0.299(11). 


I. INTRODUCTION 

Precision tests of the Standard Model (SM) seek to find discrepancies that may indicate 
the presence of new physics. The SM reqnirement of a nnitary Cabibbo-Kobayashi-Maskawa 
(CKM) weak mixing matrix provides a good opportnnity for snch a test. The nnitarity- 
triangle test checks the orthogonality of the first and third rows of the CKM matrix. It 
reqnires consistency between resnlts extracted from the experimental measurements and 
theoretical calculations of a wide variety of flavor- and CP-violating observables. Although 
recent results have been roughly consistent with unitarity [1, 2], some disagreements persist 
and require further attention. The CKM parameter \Vcb\ plays an important role in the 
unitarity triangle test, since it normalizes the lengths of sides of the triangle and contributes 
to tension in the unitarity constraint. 

The SM parameter \Vcb\ is determined through the combination of an experimental mea¬ 
surement of a branching fraction and the theoretical calculation of the underlying hadronic 
form factor(s). There are two common approaches [3] using (1) the exclusive processes 
B —)■ Dtv and B —>■ D*iV with lattice-QCD determinations of the relevant hadronic form 
factors [4, 5] or (2) the inclusive decay B —)■ XciR to any charm-containing hnal state Xc and 
the operator-product and heavy-quark expansions to describe strong-interaction effects [6]. 
It is troublesome that the most recent results for \Vcb\ from these exclusive and inclusive 
determinations disagree at the 3a level [3, 5]. It is unlikely that this difference is due to new 
physics effects [7], and further work is needed to refine the determinations. 

Reducing the error in the determination of \ Vcb\ requires both experimental and theoretical 
effort. Recent work by the BaBar collaboration [8] has provided better measurements of 
the decay rate. The latest results from the Belle collaboration for this process are still 
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preliminary [9]. Further improvements will come from Belle II. In this work we improve 
the exclusive determination of \Vcb\ from the decay process B —)■ Dtu by providing the 
hrst unquenched lattice-QCD calculation of the relevant form factors with a complete error 
budget and small statistical and systematic errors. 

Traditionally, experimental measurements are hrst extrapolated to zero recoil, where the 
comparison with theory to obtain \Vcb\ occurs, using a parameterization of the momentum 
dependence from Caprini, Lellouch, and Neubert (CLN) [10]. Indeed lattice calculations at 
zero recoil momentum are simpler, and for the exclusive process B —)■ D*iV this method 
yields a very precise determination of \Vcb\, which is described in our companion work [5]. 
However, due to the more severe phase space suppression of the B —)• Diu rate near zero 
recoil (compared with B —)■ D*tu) the extrapolation of the experimental data to zero recoil 
is a source of signihcant uncertainty. This results in determinations of |V"c6| from B —Dtu 
that are less precise than they have to be. Here we calculate the form factors for B —Diu 
for a range of recoil momenta and parameterize their dependence on momentum transfer 
using the model-independent ^-expansion of Boyd, Grinstein and Lebed [11]. We ht the 
experimental and lattice data together as a function of momentum transfer and determine 
\Vcb\ from the relative normalization over the entire range of recoil momenta. 

Where previous calculations of this process at nonzero recoil ignored effects of sea 
quarks [12], ours includes them. The use of asqtad-improved staggered fermions and 
improved Wilson (“clover”) quarks reduces lattice discretization errors. A preliminary 
determination of |I4fe| from B —)■ Dtu with a very small subset of the present asqtad en¬ 
sembles was presented in Ref. [13]. That calculation was done only at zero recoil and used 
only lattices with spacing approximately 0.12 fm, so a continuum extrapolation was not 
possible. The broad range of lattice spacings and sea-quark-mass ratios in our current study 
gives conhdence in the extrapolation to physical quark masses and zero lattice spacing. 
More recently, in a related project of ours [14], the ratio of form factors for Bg —?■ Dgli' to 
B —)■ Diu decays was obtained using only four asqtad ensembles, i.e., with two different 
light sea-quark masses at two lattice spacings. This data set was also used to obtain the 
hrst Standard-Model prediction for R{D) = B{B —)■ Dti')/B{B — D^v) from unquenched 
lattice QCD in Ref. [15]. The present work uses all fourteen ensembles and uses a slightly 
different analysis. These are the hrst such calculations to combine all of the ingredients 
listed above. 

Preliminary results for the B ^ D form factors from this project were presented in [4]. 
The hnal analysis presented here includes a more sophisticated treatment of the matching 
factors as well as more rehned estimates for the renormalization and heavy-quark discretiza¬ 
tion errors. 

This article is organized as follows. In Sec. H we review the formalism and our strategy 
for extracting the form factors at nonzero recoil. In Sec. HI we give details of the ensembles 
and simulation and discuss our determination of the form factors and the chiral-continuum 
extrapolation. We discuss systematic errors in Sec. IV. In Sec. V we present hts to our 
lattice data for the two form factors /_|_ and /o and a joint ht that combines our lattice data 
with the 2009 BaBar measurements [8], leading, hnally, to our result for |Wb|- We discuss 
our results in Sec. VI. Appendix A discusses technical details regarding the tuning of the 
bare-quark masses. Appendix B derives the pattern of heavy-quark discretization ehects 
and discusses some details of matching lattice gauge theory with heavy quarks to continuum 
QCD. 
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II. FORM FACTORS 


A. Continuum form factors 


The hadronic interaction in the process B —)■ D^T> is determined by the transition matrix 
element of the vector current = c'y^b, which is conventionally decomposed in terms of 
the vector and scalar form factors /+(g^) and /o(5'^) as 


{D{pn)mB{pB)) = U{q^ 


{Pb+PdT 


Ml-Ml 


qz 


( 2 . 1 ) 


Here pb and po are the momenta of the B and D mesons, Mb and Mb are the respective 
masses, and q = Pb — Pd is the momentum transferred to the leptons. In the approximation 
that the masses of the leptons i = e, p, Ve, are much smaller than the B and D mass 
difference Mb — Mb, the differential decay rate is 

T(b ^ DtU) = |r7Ewr^^d|!M(„2 _ + rmwf , (2,2) 

where |r/EwP accounts for electroweak corrections discussed below, Gp the Fermi weak 
decay constant, \Vch\ is the desired CKM matrix element, w = v ■ v' is the recoil parameter, 
V = pb/Mb and v' = pb/Mb are the hadronic velocities, and Q is related to /+ through 




(1 + r)^ 

4r 


0(wf . 


(2.3) 


for r = Mb/Mb = 0.354. 

The alternative parameterization in terms of the form factors h+ and h- is convenient in 
heavy-quark effective theory (HQET) and heavy-light meson chiral perturbation theory: 


{D{pB)mB{pB)) 


\/ MbMb 

These form factors are related to /+ and /o through 


= h+{w){v + v'Y + h-{w){v — v'Y ■ 


/o(g^) 


2v^ 


[(1 -h r)h+{w) - (1 - r)h_{w)] , 


w + 1 
1 -|- r 


h+{w) 


w — 1 
1 — r 


h_{w) 


(2.4) 


(2.5) 

( 2 . 6 ) 


where = Ml -|- Ml — 2 wMbMb- We note, also, the kinematic constraint /+(0) = /o(0) 
at = 0, which corresponds to tc = {Ml -|- MI)/{2MbMb) ~ 1.59. We also have 


G{w) = h+{w) 


1 -|- r 


h_{w) 


(2.7) 


B. Form factors from lattice matrix elements 

We use the local Fermilab-improved vector current for the quark transition x ^ y 

( 2 . 8 ) 
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where the subscripts denote flavor, \1/ is the “rotated” held [16] 

d' = (1 + di7 ■ Aat) V’, (2.9) 

and V’ is the heavy-quark held in the action. The lattice current is related to the 
continuum current through 

( 2 , 10 ) 

where “=” denotes the equality of matrix elements. Following [17, 18], we define the correc¬ 
tion matching factor as the double ratio of matching factors for flavor off-diagonal currents 
to those for flavor-diagonal currents: 


2 


Zy^ Zy^ 

^cb ^he 
yfjb 


( 2 . 11 ) 


where Pyi^ = 1 + 47rQ;s(g*)pyi, -|- 0(as(q*)^) is determined to one-loop order in lattice per¬ 
turbation theory [18]. It is found to be quite close to 1 because of cancellations in the ratio 
of similar quantities, including cancellations of tadpole diagrams. The truncation error is 
expected to be small because as(g* = 2/a) ~ 0.2. 

The matching factor pYt^{w) depends, in principle, upon the velocity transfer w. At 
present we have calculated only pv^iX) the quark masses and lattice spacings in our 
project. Calculation of the spatial correction pYi is more difficult because, even for zero 
recoil, one must calculate it for nonzero momentum. Thus we have calculated Pv^X) only 
for the simpler case rrica = 0, but our lack of knowledge of the rric dependence of the one-loop 
correction to py* makes only a small contribution to our final uncertainty. The w dependence 
of pYi is also unavailable. Below we note where these issues arise. 

To compute the form factors h._|_ and h_ at arbitrary recoil, we need the lattice matrix 
elements of both the temporal and spatial vector currents, and V. In practice, we 
use ratios of lattice correlators in which the flavor-conserving renormalization factors are 
automatically included, as discussed below. These ratios also suppress statistical fluctuations 
and systematic errors. The remaining correction factors py4 and pyi are applied after fitting 
the ratios. We apply this correction in Sec. Ill F. 

Our calculation is done in the S-meson rest frame for any recoil D-meson momentum p. 
We compute the double ratio 

O (P(0)|V;t|g(0))(g(0)lnl|C(0)) , , 

+ (D(0)|V'i|D(0)){i?(0)|V;‘|B(0)) ' ■ ' 


and the single ratios 


r Ir,) - <'D(p)l''ts(0)) 
(D(0)|I/4|5 (o)) ’ 

(2.13) 

_ {D(p)|V'|i?(0)) 

{D{p)\V^\B{0))^ 

(2.14) 

- {C(p)|V'*|D(0)> ^ 

(2.15) 


Note that Q+(p) is the ratio of B —?■ V matrix elements at nonzero and zero recoil, and 
that Xf{p) is computed only from the flavor-diagonal transition D ^ D. As spelled out 


5 







below, we use R+{p), Q+{p), and R~{p) to obtain h^{w) and h_{w), and Xf{p) to obtain 
the recoil w. The flavor-conserving renormalization factors and ZyA^ cancel exactly in 
the double ratio -R+, which was introduced by Hashimoto et al. and used to obtain the 
B —> D^l> form factor at zero recoil in quenched lattice QCD [19]. 

From Eq. (B15), the 3-vector Xf yields the velocity without any matching ambiguities: 


V 


Xf = 


tn -I- 1 


Because w'^ = 1 + (when the initial meson is at rest), one finds 


(2.16) 


w{p) 


1 - x){p) ■ 


(2.17) 


Thus, even the kinematic variable w is determined dynamically from a ratio of matrix 
elements. 

The other ratios require matching factors. It is convenient to define 


R+ = pI.(1)R+ , (2.18) 

C+(p) = ((d^Q+(p), (2.19) 

TC-(p) = ^4fl-(p). (2.20) 

Pv*[w) 


We derive these factors and discuss how we handle them in Appendix B. Note that 77.+ 
reduces to 

\/TI+ = h+(l) -|- matching & discretization errors . (2.21) 

Also, Q+(0) = Q+(0) = 1 by construction. We then can obtain h+ and h- from 


h+ (w(p)) = v^Q+(p) 
h- {w{p)) = ^/n^Q+{p) 


l-iz^{p) -Xfip)^ , 

_ 7l-{p)-Xf{p) 

xj{p) 


( 2 . 22 ) 

(2.23) 


as in Eq. (2.21) up to matching and discretization errors. 


III. ANALYSIS 

A. Lattice action and parameters 

Our calculation uses fourteen ensembles of gauge-held configurations generated in the 
presence of 2-1-1 flavors of asqtad-improved staggered sea quarks by the MILC collaboration 
[20]. Ensembles are indicated graphically in Fig. 1, and they are tabulated in Table I. There 
are four lattice spacings, a ~ 0.12 fm, 0.09 fm, 0.06 fm, and 0.045 fm, and light sea-quark 
to strange sea-quark mass ratios m'/m'g ranging from 0.4 to 0.05. The strange sea-quark 
mass is set approximately to its physical value. For the light valence quarks we use the 
asqtad action. Light-quark propagators are converted to improved “naive” propagators as 
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FIG. 1. (color online) Range of lattice spacings and light-quark masses used here. The area of 
each disk is proportional to the number of conhgurations in the ensemble. 


in Ref. [21] to implement the standard Dirac spin algebra. In this study, masses of the light 
valence quarks are always equal to the sea-quark masses. For the heavy valence quarks we 
use the Fermilab interpretation of the clover action with the parameters listed in Table II. 

Two-point and three-point correlators are computed from four equally-spaced source 
times per conhguration, but with random offsets in time and space to reduce correlations 
between successive gauge-held conhgurations within an ensemble. We performed a blocking 
study to look for residual autocorrelations, and found that the statistical errors did not 
change signihcantly with block size. Thus we do not block the data in this work. The 
masses of the heavy valence quarks were tuned so that the kinetic masses of the Dg and 
Bg mesons were equal to their physical values. A detailed discussion of tuning is given in 
the appendix of Ref. [5], where we show that we get good agreement between the lattice 
values of the Dg and Bg hyperhne splittings and their experimental values. The simulation 
values of the heavy-quark masses are not quite the same as our best-tuned values, which 
were determined a posteriori. Post-simulation adjustment for heavy-quark-mass tuning is 
described in Sec. HID. 

After hxing the lattices to Coulomb gauge, two types of interpolating operators for the 
D meson are used, namely, a local operator and a smeared operator based on a Richardson 
IS wave function [22]. For the B meson we use only the IS operator. These two operators 
have different overlap with excited states, so computing both helps us remove excited-state 
contributions. We generate three-point functions in a standard way by hxing the position of 
the D and B mesons to a separation T in imaginary time and then varying the time t of the 
vector current. Calculations at two adjacent time separations T are carried out in each case 
to control the effects of oscillating staggered-fermion propagators. We rotate the heavy-quark 
helds as in Eq. (2.9) using the tadpole-improved tree-level values for di listed in Table II, 
so that the vector current is tree-level improved. Calculations are made at several choices 
of three-momentum. In units of 27i/L, for this study we use hve momenta (0,0,0), (0,0,1), 
(0,1,1), (1,1,1), and (2,0,0). Results at larger momenta tend to have signihcantly larger 
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TABLE I. Parameters of the lattice-gauge-field ensembles. The columns from left to right are 
the approximate lattice spacing in fm, the bare sea-quark masses in lattice units am'/am'g, the 
lightest pseudoscalar in MeV, the root-mean-square (RMS) mass of the pion taste multiplet in 
MeV, the dimensionless factor M^L, the dimensions of the lattice in lattice units, the number 
of configurations in each ensemble (four sources each), and the tadpole-improvement factor uq 
(obtained from the average plaquette). 


a (fm) 

am'/am'g 

(MeV) 

j^RMS (MeV) 

M^L 

Lattice size 

Configs 

Uo 

Ri 0.12 

0.02/0.05 

560 

670 

6.2 

20^ X 64 

2052 

0.8688 


0.01/0.05 

390 

540 

4.5 

20^ X 64 

2259 

0.8677 


0.007/0.05 

320 

500 

3.8 

20^ X 64 

2110 

0.8678 


0.005/0.05 

270 

470 

3.8 

243 X 64 

2099 

0.8678 

« 0.09 

0.0124/0.031 

500 

550 

5.8 

28^ X 96 

1996 

0.8788 


0.0062/0.031 

350 

420 

4.1 

28^ X 96 

1931 

0.8782 


0.00465/0.031 

310 

380 

4.1 

32^ X 96 

984 

0.8781 


0.0031/0.031 

250 

330 

4.2 

40^ X 96 

1015 

0.8779 


0.00155/0.031 

180 

280 

4.8 

64^ X 96 

791 

0.877805 

Ri 0.06 

0.0072/0.018 

450 

470 

6.3 

48^ X 144 

593 

0.8881 


0.0036/0.018 

320 

340 

4.5 

48^ X 144 

673 

0.88788 


0.0025/0.018 

260 

290 

4.4 

56^ X 144 

801 

0.88776 


0.0018/0.018 

220 

260 

4.3 

64^ X 144 

827 

0.88764 

« 0.045 

0.0028/0.014 

320 

330 

4.6 

64^ X 192 

801 

0.89511 


statistical errors, and also suffer from larger momentum-dependent discretization errors. In 
the two-point correlator these momenta are projected at the sink and in the three-point 
correlator, at the current. In the latter case the three-momentum of the B meson is set to 
zero. 


B. Fitting strategy 

We need both two-point and three-point correlation functions to construct the form factor 
introduced in Sec. II B. We use interpolating operators Oxa{p,t) of spatial momentum p 
and time t with X G {B,D} and a G {15, d}. The notation d signifies a delta function 
(point) source, while 15 denotes a IS Richardson wavefunction. See Ref. [23] for details. 
The correlation functions can be expressed in terms of operator matrix elements: 

t) = (oL(p, 0)OxJp, t)) , (3.1) 

= (OP(-p,0)V'»‘(p,()Ox«(0,T)) , (3.2) 

where T is the imaginary time separation between the B and D mesons. 

The spectral decomposition of the two-point correlator is 

C.2p.,x««,p_ ^ ii,.(t) v'Zx.MZp M + eM-En{pm - «))l . 

(3.3) 
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TABLE II. Parameters of the heavy valence quarks. The approximate lattice spacing and bare 
sea-quark masses in the first two columns identify the ensemble. The remaining columns show the 
coefficient of the clover term in the SW action cswj the bare hopping-parameter k, and the rotation 
parameter in the current di. The primes on k distinguish the simulation from the physical values. 


PS a (fm) 

am!/am'g 

csw 


dib 


die 

0.12 

0.02/0.05 

1.525 

0.0918 

0.09439 

0.1259 

0.07539 

0.12 

0.01/0.05 

1.531 

0.0901 

0.09334 

0.1254 

0.07724 

0.12 

0.007/0.05 

1.530 

0.0901 

0.09332 

0.1254 

0.07731 

0.12 

0.005/0.05 

1.530 

0.0901 

0.09332 

0.1254 

0.07733 

0.09 

0.0124/0.031 

1.473 

0.0982 

0.09681 

0.1277 

0.06420 

0.09 

0.0062/0.031 

1.476 

0.0979 

0.09677 

0.1276 

0.06482 

0.09 

0.00465/0.031 

1.477 

0.0977 

0.09671 

0.1275 

0.06523 

0.09 

0.0031/0.031 

1.478 

0.0976 

0.09669 

0.1275 

0.06537 

0.09 

0.00155/0.031 

1.4784 

0.0976 

0.09669 

0.1275 

0.06543 

0.06 

0.0072/0.018 

1.4276 

0.1048 

0.09636 

0.1295 

0.05078 

0.06 

0.0036/0.018 

1.4287 

0.1052 

0.09631 

0.1296 

0.05055 

0.06 

0.0025/0.018 

1.4293 

0.1052 

0.09633 

0.1296 

0.05070 

0.06 

0.0018/0.018 

1.4298 

0.1052 

0.09635 

0.1296 

0.05076 

0.045 

0.0028/0.014 

1.3943 

0.1143 

0.08864 

0.1310 

0.03842 


where there are either nonoscillating terms with Sn{t) = 1 or staggered-fermion opposite- 
parity oscillating terms Snif) = —(—!)*, iV* is the lattice extent in time, and Zxa,n is the 
overlap coefficient. For the three-point function, the decomposition is similar: 


^-En{p)t 


n,m V^KiP) 


(y6,n(p)|I/^|Xa,m(0))(3.4) 


X 


where we have assumed t < T -C At, so we may neglect wraparound terms with t ^ Nt — t 
and T — t ^ Nt — {T — t). 

The double ratio R+ can be calculated very precisely from 


R+,b(t, T) 


^3pt,D6— 


(M. 

(0,f)‘ 


(3.5) 


This quantity depends on t, T and the A-meson interpolating operator, labeled by b. The 
dependence arises from contributions from excited states and opposite-parity oscillating 
states. As in Refs. [5, 24, 25] we suppress contributions from oscillating states by averaging 


R+,b(ty 


'^R+,b{t^T) + + 1) + -k 1,T -k 1) . 


(3.6) 


We drop the bar henceforth. We use a similar method for the other three-point correlation 
functions. We find that the suppression of oscillating states for B ^ D correlators is similar 
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to that of our previous work on i? —)■ D*ii> [5]. In particular, the contribution from the hrst 
oscillating B- and D-meson excited states, which does not itself oscillate in time, is reduced 
by a factor of ~ 5-12 using the average in Eq. (3.6), where greater suppression occurs for 
hner lattice spacings. 

For large t and T — t, excited-state contributions are negligible, giving the desired result. 


R+,b(t, T) —)■ R 


(3.7) 


as a plateau in the ratio vs. t, as illustrated in Fig. 2. The leading corrections to the 
plateau arise from contributions from the hrst excited D- and S-meson states. For large 
t and T — t, their contributions to the correlator double ratio fall off as exp[—AMf] and 
exp[—AM(T — f)], where AM = AMb or AMo, the splitting between the ground state 
and hrst excited state of the B- and H-mesons, respectively. Since they are both small, for 
htting the ratio, we use the approximation 

R+^b{t,T) ^ R+ + Ai?^,bexp(-AMM) + 5R^,bexp[-AMB(T - t)] 

+ C'ij_^,feexp(-AMsf) DR^^bexp[-AMD(T - f)] exp^AMst ), (3.8) 

However, since AMb ~ AMb we construct the ht model from only the i?+. A, and B terms. 
Similarly, we introduce a time- and interpolating-operator-dependent ratio 


Q+ApA,t) 


gf t) EoZm(0) 

MoZnip) 


(3.9) 


In this ratio the plateau takes on the value Q+{p) introduced in Eq. (2.13). Again, the 
leading corrections to the plateau arise from contributions from the hrst excited D- and 
H-meson states. For large t and T — t, their contributions to the correlator ratio fall oh as 
exp[—AF/jf] and exp[—AM 5 (T — f)]. Where they are both small, for htting the ratio, we 
use the approximation 

(5+,fe(p,f,T) ^ Q+{p) exp{6mt) -F AQ_^,fe(p) exp(-AFM) 

-h Bq^^bip) exp[-AMDt] Cg^fiip) exp[-AMB(T - t)]. (3.10) 

The parameter 6m vanishes when the exponential factor in Eq. (3.9) cancels the time depen¬ 
dence in the three-point functions, as it should. Since there may be slight diherences in the 
determination of the masses from the three-point and two-point functions, the cancellation 
might not be perfect. Therefore, we introduce 5m as a constrained htting parameter. The 
prior constraint is centered at zero and it has a width determined from the small statistical 
error in the two-point-htted energies. In practice, the values of aSm, are typically of order 

10-h 

For Rl_ fj{p, t, T) we form the ratio 


Rt,bipA,T) 




(3.11) 


SO that for large t and T — t we have R_^b{p,t,T) —)■ R_{p). 
use the ratio 

^3pt,Db^Dbf 

A (b f T) = — _ 


Similarly, for xf{p,t,T), we 

(3.12) 
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TABLE III. Comparison of ground-state energies Ed, excited-state energies E'^, and ground-state 
overlap factors and Zis for the (2 -|- 2)-state and (3 -|- 3)-state two-point-correlator fits for the 
D meson on the a ~ 0.12 fm, m' = O.ldm^ ensemble. In all cases the (2 -|- 2)-state fitting range is 
[6, 16] and the (3 -|- 3)-state htting range is [4, 23]. We use (3-|-3)-state hts for the analysis; the 
(2-|-2)-state fits just provide a check of systematic effects. 


P (2vr/L) 

Ground state {oEd) 
2+2 3+3 

1®* excited state {aEfi) 
2+2 3+3 

X^/df 

Zis,is Zd,d 2 + 2 3 + 3 

000 

0.9566(6) 0.9566(7) 

1.54(3) 

1.41(4) 

4.045(29) 0.0785(7) 18.9/18 35.5/37 

100 

1.0013(10) 1.0017(9) 

1.55(2) 

1.39(4) 

2.912(36) 0.0741(10) 16.5/18 46.6/37 

no 

1.0436(15) 1.0433(12) 

1.56(2) 

1.41(3) 

2.149(38) 0.0704(13) 16.5/18 36.9/37 

111 

1.0838(21) 1.0831(15) 

1.60(2) 

1.45(3) 

1.628(40) 0.0673(18) 22.6/18 39.3/37 

200 

1.1206(31) 1.1172(23) 

1.60(3) 

1.48(4) 

1.279(46) 0.0658(25) 17.4/18 46.9/37 


TABLE IV. Fit ranges [tminj ^max] for two-point and three-point functions. They are chosen to be 
approximately similar in physical units and independent of sea-quark masses with one exception: 
for the case a ~ 0.12 fm and w!jm!^ = 0.1, the two-point range was [3,23]. 


Ri a (fm) 

two-point three-point 

0.12 

[4,23] 

[2,10] 

0.09 

[5,33] 

[2,15] 

0.06 

[7,45] 

[4,18] 

0.045 

[11,80] 

[7,24] 


so that Xf^b{p,t,T) —)■ Xf{p). To fit the time dependence of R-^b{p,t,T) we use 

R-^b{p,t,T) ^ R-{p) + AR_^b{p) exp(-AEi 5 t) -b exp[-AMs(T - t)], (3.13) 

and for the time dependence of Xf{p,t,T), we use the same form, except replacing AMr 
with AMr. 


C. Correlator Fitting 

We obtain the lattice form factors via a two-step procedure. First, we £t the B- and D- 
meson two-point correlators to obtain the energies and overlap factors. Then we use these 
results as constraints with Bayesian priors in the three-point hts. Errors in the resulting form 
factors h+ and are determined from a complete single-elimination jackknife procedure. 


1. Two-point correlator fits 

The two-point functions are constructed from both a local and a smeared interpolating 
operator. They are £t simultaneously to determine the ground- and excited-state energies. 
We include oscillating and nonoscillating states in pairs and test the stability of the hts 
by comparing results with 2-1-2 and 3-1-3 states. An example is shown in Table III for the 
a ~ 0.12 fm, m' = 0.14m( ensemble. For this case we choose a fit range of [4,23] with 
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FIG. 2. Sample joint three-point function fits for determining the ratios and Q^{p) (left), 
(middle) and xj.(p) (right), for lattice momentum p = (1,0,0). Data shown are for the 
a ~ 0.06 fm, m' = 0.14m^ ensemble with B-D separation T = 24,25. Values are plotted against 
the time t of the vector-current insertion. Data points at the left and right extremities are not 
included in the fit. A color (gray-scale) change indicates which points are included in the fit. 
Black lines indicate the upper and lower Icr range of the ground-state contribution. Best fit lines 
are shown in red (gray). For Q+{p) the “plateau” is slanted because of the factor exp(<5mt) in 
Eq. (3.10). 


3-1-3 states. Results for the energy and overlap factor for that range agree with hts in the 
range [6,16] with 2-1-2 states. For the analysis, we use hts with 3-1-3 states; the (2-1-2)- 
state hts provide a check of systematic ehects from excited-state contamination. We select 
approximately the same ht ranges in physical units for all ensembles, as shown in Table IV. 


2. Three-point correlator fits 


To determine the nonzero-recoil form factor R^{p), we ht three ratios simultaneously: 
the double ratio for the IS source from Eq. (3.5) and the local- and smeared-source ratios 
Q+^d{p,t,T) and Q+^is{p,t,T) from Eq. (3.9). Because the ht model [Eq. (3.10)] includes 
ehects of the same hrst-excited states that occur in the two-point functions, we use the 
two-point-ht values for these states to set priors for Sm, SEd, and 6Mb- The best-ht values 
are used as the central values and their errors as the widths of the Gaussian priors. Fit 
ranges are chosen for stability. We use the same range for all three-point correlators in a 
given ensemble. The ranges are listed in Table IV, and a sample three-point ht is plotted in 
Fig. 2, left. 

For Xf^b{p, t, T) and R-^b{p, t, T), we ht values for both local and smeared sources jointly 
with the htting form of Eq. (3.13). Sample three-point hts are plotted in the middle and 
right panels of Fig. 2. Then, having determined all the needed quantities, we calculate w, 
h^{w) and h_{w) from Eqs. (2.17) and (2.22)-(2.23) for each momentum p and ensemble. 
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D. Heavy-quark-mass adjustment 


We adjust the bare masses of the b and c quarks so that the kinetic masses of the Dg 
and Bg mesons obtain their physical values. When computing the two-point and three- 
point correlators, we used good estimates of these quark masses. By the end of the data 
generation, we could obtain better estimates via the procedure described in [5]. 

Because there are small differences between the simulation values and hnal, tuned values, 
an adjustment of the form factors is required. Details are given in Appendix A. To obtain 
the adjustment we computed a full set of correlation functions on one of our ensembles with 
a few heavy-quark masses close to the tuned value and use these results to calculate the 
slopes of the form factors with respect to the quark masses. These results and the known 
corrections then give the needed small adjustments tabulated in Table XII. The size of the 
heavy-quark mass corrections to h_|_ (h_) range from 0 to 0.2% (0 to 2%). Small errors arise 
both from uncertainties in the tuned quark masses and uncertainties in the determination 
of the slopes. 


E. Current renormalization 

Here we summarize the procedure for matching the lattice matrix elements to the con¬ 
tinuum. The three-point hts yield ratios in which the flavor-diagonal factors ZyA ZyA from 
Eq. (2.11) cancel. Thus, to normalize the form factors to continuum conventions, we only 
have to apply the flavor off-diagonal factors pyti. as in Eqs. (2.18)-(2.20). Matching fac¬ 
tors with two heavy quarks depend on the recoil tc, but the w dependence is not available. 
Even so, we can obtain some information by considering the limit m 2 cCb 1, where the w 
dependence goes away. For this reason, each of the matching factors in Eqs. (2.18)-(2.20) 
requires somewhat different treatment. Appendix B provides further details on the matching 
calculations. 

The calculation of the zero-recoil matching factor pyA{l) needed to renormalize R+ is 
completely analogous to that of the axial-vector matching factor used in Ref. [5] . Following 
Ref. [18], we compute it to one-loop order in perturbation theory, 

pyA{l) = l + ay{q*)p%{l), (3.14) 

where ay{q*) is the QCD coupling in the V scheme [26], evaluated here at the scale q* = 2/a. 
The result for each ensemble is listed in Table V. 

For the matching factor pyA^w)/pyA^l), we note that, by construction, the one-loop coef- 
hcient must be proportional to tc — 1. Moreover, for m 2 c(i 1, which holds on the two hnest 
lattices, one may treat the charm quark as a light quark [18], using the HQET formalism 
for heavy-light currents [27]. The w dependence goes away in this limit, so the one-loop co¬ 
efficient must also be proportional to m 2 ca. Thus, pyA^w)/pyA^l) = 1 + 0 {ag{w — l)m 2 ca), 
where the coefficient of the one-loop correction is not known. In our analysis, we take 
pyA^w)/pyA^l) = 1 and include the estimated size of the one-loop correction as a w- 
dependent uncertainty. 

For the matching factor pyi{w)/py a (w), we can take the heavy-light theory a step further 
and calculate the matching explictly for m 2 ca <C 1. The calculation does not depend on w. 
The resulting values for pyi{w)/pyA[w) as in Eq. (3.14) are listed in Table V. The error in 
the one-loop coefficient introduced by taking the limit m 2 ca —)■ 0 is proportional to agm 2 ca 
with a, presumably, mild w dependence, and is again included as an uncertainty. 
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TABLE V. One-loop estimates of the matching factors for the lattice ensembles in this study. 
Shown are the approximate lattice spacing in fm, the sea-quark mass ratio m!jm'g, the tuned k 
values of the charm and bottom quarks [5], the strong coupling in the E-scheme evaluated at 
q* = 2/a, and the zero-recoil factors Pv^{^) and pyi{w)/pyi{w) on that ensemble. The first error 
in each tuned k value is statistical, and the second reflects the uncertainty in the lattice scale 
determination [5]. The correction factors are evaluated at the tuned heavy-quark masses except 
for pyijpyi^ which is evaluated at rrica = 0. The systematic uncertainties in the p factors are 
discussed in Sec. IV and Appendix B. 


« a (fm) 

m'/m'g 

Kc 

Kb 

OiviQ* = 2/a) 

/5y4(l) 

Pvi{w)/pv4{w) 

0.12 

0.4 

0.12452(15)(16) 

0.0879(9)(3) 

0.3047 

1.025105 

0.892347 

0.12 

0.2 

0.12423(15)(16) 

0.0868(9)(3) 

0.3108 

1.026472 

0.888051 

0.12 

0.14 

0.12423(15)(16) 

0.0868(9)(3) 

0.3102 

1.026395 

0.888248 

0.12 

0.1 

0.12423(15)(16) 

0.0868(9)(3) 

0.3102 

1.026388 

0.888241 

0.09 

0.4 

0.12737(9)(14) 

0.0972(7) (3) 

0.2582 

1.015603 

0.924664 

0.09 

0.2 

0.12722(9)(14) 

0.0967(7) (3) 

0.2607 

1.016080 

0.923051 

0.09 

0.15 

0.12718(9)(14) 

0.0966(7) (3) 

0.2611 

1.016160 

0.922757 

0.09 

0.1 

0.12714(9)(14) 

0.0965(7)(3) 

0.2619 

1.016259 

0.922319 

0.09 

0.05 

0.12710(9)(14) 

0.0964(7) (3) 

0.2623 

1.016340 

0.922022 

0.06 

0.4 

0.12964(4)(11) 

0.1054(5)(2) 

0.2238 

1.008792 

0.947870 

0.06 

0.2 

0.12960(4)(11) 

0.1052(5)(2) 

0.2245 

1.008945 

0.947361 

0.06 

0.14 

0.12957(4)(11) 

0.1051(5)(2) 

0.2249 

1.009017 

0.947085 

0.06 

0.1 

0.12955(4)(11) 

0.1050(5)(2) 

0.2253 

1.009098 

0.946829 

0.045 

0.2 

0.130921(16)(70) 

0.1116(3)(2) 

0.2013 

1.004566 

0.962520 




FIG. 3. Global fit of all data for the form factors (left) and h- (right) vs. recoil w. The blue 
(shaded) band gives the la confidence range for the continuum extrapolation at physical quark 
masses. Fit errors include statistics, matching and truncation of the chiral expansion. The legend 
in the left figure gives the color convention for the lattice spacing, and, in the right, it gives the 
shape convention for the sea-quark mass ratio. 
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F. Chiral-continuum extrapolation 


The resulting form factors h_|_ and h_, after applying the k corrections and renormalization 
factors, are shown in Fig. 3. As can be seen, the dependence of h_|_ on lattice spacing a and 
light-quark-mass ratio m'/m'^ is quite mild. The form factors must be extrapolated to the 
physical average value of the up and down quark mass rim and zero lattice spacing (a —)• 0) 
(the physical point). 

To this end we fit both form factors to the following expressions: 


/z.+ (a, m', m'g, tc) = 1 + 


+ Ca 


„2 1 ^ I 1^2 


mi 


- P+iw - 1 ) + k+{w - 1 )^ Ci^+Xi + Ca,+Xa2 + Ca,a,+Xa2 


+xiXa^ + C 2 ,+x^ + w, A^) (3.15) 


h_(a, m', m' w) = - p‘^{w — 1) + k_{w — 1)^ + Ci _xi + _a;a2 

me 

-f Ca,a-xl2 + Ca,m-XlXa2 + C2-xj , (3.16) 

which contain the correct dependence on the light and strange-quark masses, lattice spac¬ 
ing, and recoil w at next-to-leading order (NLO) in chiral perturbation theory. The chiral 
logarithm term, denoted “logSsu( 3 )”, contains non-analytic dependence upon the pion and 
kaon masses (or equivalently m' and m(). It comes from a staggered-fermion version of the 
one-loop continuum result of Chow and Wise [28] that includes taste-breaking discretiza¬ 
tion effects [29]. The explicit expression for logSsu( 3 ) is given in the Appendix of Ref. [14]. 
The coefficient of the logarithm term is predicted in yPT, but depends upon the value of 
the D*-D-tt coupling, gD*DTT, which is not known precisely. We allow gD*D-K to vary in the 
fit, but constrain it with a Gaussian prior 0.53 ± 0.08, motivated by the spread of experi¬ 
mental [30-32] and recent lattice-QCD results [33-38]. The analytic terms depend on the 
light spectator-quark mass through xi = 2BQm!/fl) and on the lattice spacing through 
Xq? = [o/(47r/^ri)]^, which, according to xPT power counting, are expected to have co¬ 
efficients of order 1 [23]. The NLO expression is supplemented by next-to-next-to-leading 
order (NNLO) analytic terms in the light-quark mass m' and lattice spacing to incorporate 
the error from the truncation of the chiral expansion, as explained below, and by terms 
analytic in (w — 1) to allow interpolation in w at nonzero recoil. We do not include analytic 
functions of the strange sea-quark mass because (1) we do not have sufficiently varied values 
of m'g to be able to resolve any strange-quark mass dependence, (2) from yPT we expect 
the sea-quark mass dependence of the form factors to be signihcantly smaller than the light 
spectator-quark mass dependence, and (3) as discussed in Sec. IV G, we do not observe any 
strange sea-quark mass dependence within our current statistical precision. 

The statistical errors and correlations from the two-point and three-point ratio fits are 
propagated to the chiral fits using a single-elimination jackknife procedure. The strongest 
correlations are between the data for h+{w) (or h-{w)) at different w values on the same 
ensemble. The data for h^{w) and h-{w) on the same ensemble are only weakly correlated. 
Results from different ensembles are statistically independent. The fits to Eqs. (3.15) and 
(3.16) are done taking fully into account all statistical correlations. 

To test the applicability of NLO chiral perturbation theory to our data, we first fit without 
the analytic NNLO terms. The p value, p = 0.93 of the joint, exclusively NLO ht to and 
is satisfactory.^ Next we include the analytic NNLO terms with priors 0 ± 2 based on 


With Gaussian priors our p value is determined from the augmented We count degrees of freedom as 
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FIG. 4. Comparison of NLO (hatched) and NNLO (solid) chiral-continuum fits for /i_|_ (left) and 
h- (right) vs. recoil w. 

expectations from xPT power-counting in order to test for the effect of truncating the chiral 
perturbation series. The p value decreases slightly to 0.87. Including these terms increases 
the standard deviation by at most ~10% for h_|_ and ~30% for h_, and shifts the central 
values by much less than the hnal standard deviation, as shown in Fig. 4. The statistical 
errors then can be safely assumed to include the systematic error of the truncation. We 
therefore use this £t including NNLO terms to obtain our preferred value for the form 
factors at the physical point. The results of the extrapolation with propagated statistical 
errors are shown as bands in Fig. 3. 

In heavy-quark effective theory, Luke’s theorem states that h+{w = 1) has leading cor¬ 
rections only at second order in the inverse heavy-quark masses, namely 1/m^ and 1/ml, 
whereas has corrections at hrst order. Appendix B and Ref. [39] show how Luke’s theorem 
applies in lattice gauge theory and, hence, that one expects h- to have larger heavy-quark 
discretization errors than Indeed, we see that does have a stronger dependence on 
lattice spacing than h_|_. For the determination /+, and therefore iWftl, fhe contribution of 
h_ over the entire kinematic range is small, so the larger errors in h_ do not increase the 
overall error much. These trends in lattice spacing with fourteen ensembles are consistent 
with our previous findings with four ensembles [14]. 

We build the systematic errors from k tuning and from the matching factors into the 
chiral-continuum extrapolation by forming the combined covariance matrix for the data as 
follows: 

G,, = Clf + (3.17) 

where the hrst term is the statistical covariance, and the index i runs over all data (ensembles, 
momenta, and and h_). We denote by and the shift on the fth datum due to 
the matching and K-tuning errors, respectively. Equation (3.17) conservatively assumes 
that the matching-factor errors (or K-tuning errors) are 100% correlated between all data 
points. For the systematic errors due to our matching procedure, we have estimates for the 

the number of data points minus adjustable parameters plus the number of theoretically-motivated priors. 

Very loose priors that have no impact but to stabilize the fits are not counted. 
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FIG. 5. The form factors /+ and /o as a function of the recoil w resulting from the chiral-continuum 
fit in this study (cyan band), compared with the results from [14] (cross-hatched band). The width 
of each band indicates the la error from the chiral-continuum fit, but uncertainties from the 
matching factors are included only in the cyan bands. See the text for additional details. 


uncertainty in ^^^(l), Pv^{w)/pvi{l) and pvi{w)/pv^{w) in Eqs. (B31), (B29), and (B37), 
respectively. The form factors h±{w) change the most when pv^{w) and pvi{w) / pyi{w) are 
simultaneously shifted in opposite directions. We take the average of these two shifts as an 
estimate of the p-factor error for all h±{w) on all ensembles. For the K-tuning error, we take 
the same approach in principle, propagating the uncertainties of the intercepts and slopes 
in Appendix A to shifts 6^^'^ of the form-factor data. However, we hnd that the resulting 
are negligibly small, and we therefore set them to zero in Eq. (3.17). 

Given the chiral-continuum £t results for and h_, we construct the vector and scalar 
form factors /+ and /o using Eqs. (2.5) and (2.6). Figure 5 compares our new B —)■ Dtv form- 
factor results with those from our earlier work [14] in the w range where we have simulation 
data. The curves shown are output from the chiral-continuum extrapolation, and therefore 
include the uncertainties from statistics, the chiral-continuum extrapolation, and matching 
(for the current work); they do not include the remaining systematic uncertainties, which 
we add in quadrature a posteriori in both works. We expect that the two results are largely 
independent because they have only a small subset of overlapping data (the earlier work 
included only four ensembles), and the new work includes NNLO analytic terms in the yPT 
fit function. The results are consistent for both form factors over almost all simulated w 
values, and diverge only slightly for /q for w > 1.13. The central values of the new form 
factors are slightly higher than in [14] , primarily due to explicit inclusion of the perturbative 
correction factors pv^{w)/pv^iw) which have a bigger effect on the form factor /o than on /+. 
The total errors on the form factors in this work are similar in size to those in Ref. [14], but 
the additional ensembles used in this work enable a more detailed and reliable systematic 
error analysis as described in Sec. IV. (Reference [14] focused on form-factor ratios in which 
most of the systematic errors are suppressed.) 
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FIG. 6. Error budgets for /+ and /o as a function of the recoil w. The colored bands show the 
error contribution of each uncertainty source to the quadrature sum. The corresponding error is 
provided on the right y-axis. Our lattice simulation results are for w G [0,1.16], i.e., to the left of 
the vertical line. 


IV. SYSTEMATIC ERRORS 

In this section we discuss the sources of systematic error in the lattice determinations 
of h_|_ and and their propagation to the form factors /+ and /q. As can be seen from 
Fig. 4, the magnitude of is about 5% of for the entire range of simulated w values. 
Further, the contribution of h._ to the vector form factor /_|_ is suppressed relative to the 
contribution from by the factor (1 — r)/(l + r) = 0.477, while the contribution of h_ to 
the scalar form factor /o is exactly zero at tc = 1 and grows linearly with recoil as (w — 1). 
Thus even large percentage systematic errors in /i_ lead to only small uncertainties in /_|_ 
and /q. Figure 6 shows the momentum-dependence of the error contributions to f+{w) and 
fo{w), while Table VI provides numerical values for a representative recoil w = 1.16. 


A. Overview of systematic errors in /+ and /o 

As can be seen from Fig. 6, the dominant uncertainty in both form factors arises from 
the chiral-continuum fit, which includes contributions from statistics, matching factors, and 
higher-order terms in the chiral expansion. Although we cannot strictly disentangle the 
contributions to the error from these sources, we can estimate their sizes by repeating the 
chiral-continuum fit omitting either the errors in the matching factors or the NNLO terms 
in the chiral expansion, and take the quadrature difference of the resulting error estimates. 
The contribution from “statistics” is defined to be the error in the NLO chiral-continuum 
£t to data with no matching-factor uncertainties included. This imprecise scheme does not 
guarantee that the individual errors sum to the total fit error, but, roughly speaking, we hnd 
that the statistics, matching, and truncation uncertainties in the chiral-continuum expansion 
contribute approximately equally to the error in the full NNLO fit. Despite our incomplete 
knowledge of the matching factors, we find their contributions to the uncertainty in /+ and 
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TABLE VI. Error budget (in percent) for /+ and fo at w = 1.16, which is the largest recoil value 
used in our momentum extrapolation to the full kinematic range and determination of \Vcb\ (see 
Sec. V). The first row includes the combined error from statistics, matching, and the error from 
truncating the chiral expansion resulting from the chiral-continuum fit: errors in parentheses are 
approximate sub-parts estimated as described in the text. The total error is obtained by adding the 
individual errors in quadrature. Not explicitly shown because they are negligible are finite-volume 
effects, isospin-breaking effects, and light-quark mass tuning. 


Source 

/+(%) 

/o(%) 

Statistics-|-matching-|-xPT cont. extrap. 

1.2 

1.1 

(Statistics) 

(0.7) 

(0.7) 

(Matching) 

(0.7) 

(0.7) 

(yPT/cont. extrap.) 

(0.6) 

(0.5) 

Heavy-quark discretization 

0.4 

0.4 

Lattice scale ri 

0.2 

0.2 

Total error 

1.2 

1.1 


/o to be modest. The errors from the chiral-continuum £t are under good control for the 
range of simulated lattice recoil values, but grow rapidly for w > 1.16 where we do not have 
data. 

We add the remaining systematic uncertainties a posteriori to the chiral-continuum fit 
error. We estimate the individual contributions to the form-factor error budget in the follow¬ 
ing subsections, discussing each source in a separate subsection for clarity. In practice, only 
the heavy-quark discretization errors (Sec. IV D) and lattice-scale uncertainty (Sec. IV E) 
turn out to be signihcant. 

We assume that systematic uncertainties from heavy-quark discretization effects and 
the lattice-scale uncertainty are uncorrelated, and therefore add them in quadrature. We 
then propagate them to /+ and /o according to the linear transformation Eqs. (2.5) and 
(2.6), which depends on the recoil w, taking them to be 100% correlated between w values 
and between h+ and h_. Both the lattice-scale and heavy-quark discretization errors are 
substantially smaller than the chiral-continuum £t error, and increase only slowly with w. 


B. Matching 

The p factors in Eq. (2.11) enter in the renormalization of the components of the transition 
vector current As explained in Sec. HIE these factors are estimated in one-loop lattice 
perturbation theory to the extent that such calculations are available. As discussed near the 
end of Sec. IIIF, we build the uncertainty estimates of Eqs. (B31), (B32) and (B37) into 
the chiral-continuum £t via Eq. (3.17). 

A noteworthy feature of Table VI is the size of the matching error after the chiral- 
continuum £t. Had we omitted the errors in Eqs. (B31), (B32), and (B37) from the fitting 
function, we would have to add them a posteriori, as we did for B ^ D* at zero recoil [5]. 
Following the procedure used in Ref. [5], we would assign errors of 1.4% and 1.1% for 
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FIG. 7. The form factor h^{w) at three representative values of the recoil w as a function of the 
squared lattice spacing (a/ri)^, for all ensembles with m' = 0.2m'^. Black squares denote data 
points interpolated to the same recoil value, while orange circles denote fit values interpolated 
further, so that the light-quark masses m' and m'^ correspond to the values on the lattice with 
a ~ 0.09 fm; the orange diamonds denote the continuum limit in this case. The solid curves show 
the a dependence predicted by the HQET description of cutoff effects, with A = 450 MeV. These 
trends are shown as deviations from the a ~ 0.09 fm lattice. For details, see the discussion of 
Tables XIII and XIV in Appendix B. Note that the data and ht points reflect discretization errors 
from light quarks and gluons, as well as those from the heavy quarks. 

and /o, respectively, at w = 1.16, based on the second-finest lattice with a ~ 0.06 fm and 
its value of a* = 0.225. Incorporating the matching errors into the chiral-continuum £t, 
however, allows them to vary with lattice spacing and to be informed by the data. It is 
reasonable that the additional information reduces the uncertainty to about 0.7% for both 
/+ and fo at w = 1.16, as shown in Table VI. 


C. Light-quark and gluon discretization errors 

Our improved actions have light-quark and gluon discretization errors of order and 
[5]. As discussed in Sec. IIIF, they are already included in the ht model of Eqs. (3.15) 
and (3.16). From Table VI, the errors due to the truncation of the chiral expansion and 
extrapolation to the continuum limit are about 0.6% and 0.5% for /+ and fo, respectively, at 
w = 1.16. Using simple power-counting, we would conservatively estimate the size of generic 
light-quark and gluon discretization errors on the a ~ 0.06 fm lattice to be about 1%. The 
data for which give the dominant contribution to /+ and %, do not display signihcant 
lattice-spacing dependence. Therefore, allowing the data to constrain the possible size of 
light-quark and gluon discretization effects reduces the error. 


D. Heavy-quark discretization errors 

An important uncertainty comes from discretization errors in the lattice treatment of the 
heavy quarks. Applying the theory of heavy-quark cutoff effects developed in Refs. [18, 39], 
we estimate the size of these errors in Appendix B, providing in Tables XIII and XIV 
numerical results for the errors on h± (w) from mismatches in the lattice action and currents 


20 













for a heavy-quark scale of A = 450 MeV. The value is the same as that used in Ref. [5] , and 
here we explain why this choice is reasonable in this case too. 

In Fig. 7, we show the observed lattice-spacing dependence of our simulation data for 
h+(tc) at three recoil values on the rh' = 0.2m'^ ensembles. The raw data (black squares) are 
adjusted slightly to obtain the same w values for all a using a chiral-continuum £t with the 
p-factor and K-tuning errors turned off. Thus, the error bars shown here are statistical only. 
We also use this £t to adjust the light-quark masses to those on the m' = 0.2mA a ~ 0.09 fm 
ensemble (orange circles). (In practice, shifting the strange sea-quark mass has little impact 
on the fit points.) To compare the trend with the expected heavy-quark discretization error, 
we draw the size of the effect — dehned as the difference from a ~ 0.09 fm — predicted in 
Appendix B for A = 450 MeV. 

For all values of re. Fig. 7 shows that this estimate captures most of the discretization 
effect, given the statistical scatter. Note that the £t-interpolated (orange) points make clear 
that the trend is predominantly linear in a?. This dependence is characteristic of generic 
discretization effects of the light quarks and gluons, which are already included in the chiral- 
continuum £t model, Eqs. (3.15) and (3.16). Moreover, the heavy-quark discretization effects 
turn out to be nearly linear in a^, so they, too, are mostly absorbed by the £t model. It 
does not make sense to count this well-modeled a dependence twice by, say, inflating A to 
encompass all of the variation seen in Fig. 7. That said, we do not have an argument to 
reduce the value of A used in Tables XIII and XIV below 450 MeV. Following Ref. [5] , we base 
our hnal estimate on our next-to-smallest lattice spacing, a ~ 0.06 fm, leading to the error 
estimates in Table VI. The heavy-quark discretization error is found to be small compared 
with the chiral-continuum extrapolation error. For h+{w) it ranges from approximately 
0.15% at w = 1 to 0.35% at our largest w values. For h-{w) it is approximately 20%. 


E. Lattice-scale error 

We use the distance scale ri and the relative lattice spacing a/ri to determine the lattice 
scale a. The ratio a/ri for the ensembles in this study is known quite precisely from a £t 
to a wide range of data for the heavy-quark potential [20] . For this study we use the values 
of Ti/a presented in Table III of [5]. The continuum, physical quark-mass value of ri is 
determined from studies of the light pseudoscalar-meson spectrum and decay constants. For 
this study we use ri = 0.3117(22) fm, based on the PDG value of [23]. 

Because the form factors are dimensionless, the lattice scale enters only weakly into their 
determination via: (1) tuning the heavy-quark masses, (2) setting light-meson masses in 
the chiral logarithms, and (3) hxing the location of the continuum limit. To determine the 
error due to uncertainties in ri we see how much our results shift when we change ri by one 
standard deviation. We hnd that the changes in the form factors are smaller than 0.2%. 


F. Finite-volume corrections 

The hnite-volume effects can be estimated within NLO heavy-light meson yPT by replac¬ 
ing the loop integrals with discrete sums. The corrections to the integrals in the formulas 
appearing in B ^ D decays at zero recoil were worked out by Arndt and Lin [40]. At the 
values of quark masses and volumes at zero recoil where we have data, the effects predicted 
by yPT are less than one part in 10^. This is not a result of cancellation, but is due to 
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the fact that the chiral logarithms make only a very small contribution to the form factor. 
We did not calculate the hnite-volume corrections at nonzero recoil because the integrals 
appearing in those formulas are much more complicated, but there is no reason to expect 
these effects to be signihcantly enhanced away from the zero-recoil point. Thus, hnite-size 
effects are expected to be negligible compared with our other errors, and we do not assign 
any additional error due to them. 


G. Light-quark-mass tuning 

We extrapolate the form factors to the physical average of the up- and down-quark masses 
rim = 0.003612(126), determined from an analysis of the light pseudoscalar-meson spectrum 
and decay constants on the same ensembles [41]. Varying rim by plus and minus la in our 
chiral-continuum £t leads to relative changes of order 10“® for both form factors in the range 
of simulated recoil values. 

On some ensembles the strange sea-quark mass deviates by as much as 30% from its 
physical value. From heavy-light meson yPT, we expect the B ^ D form factors to be 
largely insensitive to sea-quark masses. Nevertheless we study the impact of the strange sea- 
quark mass by calculating the ratios in Eqs. (2.12)~(2.15) on an a ~ 0.12 fm ensemble with 
an unphysically-light strange sea quark, am'/am(, = 0.005/0.005. We do not observe any 
statistically-signihcant differences in these ratios from those on the am'/am'^ = 0.005/0.05 
ensemble. We therefore conclude that errors from mistuning the strange sea-quark mass are 
negligible within our current precision. 


H. Heavy-quark-mass tuning 

As described in Sec. Ill D, we adjust the simulation data before the chiral-continuum £t to 
account for the slight difference between the simulated bottom and charm k values and the 
physical ones, using the corrections estimated in Appendix A. The size of these corrections 
is quite small, ranging from 0 to 0.2% for and from 0 to 2% for h_. Repeating the chiral- 
continuum £t omitting the k corrections does not appreciably change the chiral-continuum 
fit result. Thus we conclude that the uncertainty in the form factors due to errors in the 
heavy-quark masses is negligible. 


I. Isospin correction 

In our calculation we have assumed that the up- and down-quark masses are equal, 
although in nature, they are not. Therefore, if we distinguish between them in calculating 
the value of the form factors at the physical point, we get a slightly different result. To 
estimate the sensitivity of /+ and % to isospin splitting of the light-quark masses, we use 
our best-£t parameters in the chiral-continuum model, and evaluate the £t function at 
the physical values of rim„ = 0.002236 and rim^ = 0.004988, instead of rim given above. 
These values are obtained by combining rim obtained on the asqtad ensembles with the ratio 
’rriu/md = 0.4482 (l^o?) obtained from the MILC Collaboration’s study of electromagnetic 
effects on the pion and kaon mass-splittings on the (2 -|- 1 -|- l)-flavor HISQ ensembles [42]. 
The relative shifts in both form factors for all simulated recoil values are of order 10“^, and 
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TABLE VII. Selected values of the form factors f+{w) and fo{w) at the physical point (syn¬ 
thetic data) and their correlations. Errors shown include statistics and all systematics added in 
quadrature. 



value 

Correlation matrix 

/+(1) /+(1.08) /+(1.16) /o(l) /o(1.08) /o(1.16) 

/+(!) 

1.1994(095) 

1.0000 0.9674 0.8812 0.8290 0.8533 0.8032 

/+(1.08) 

1.0941(104) 

1.0000 0.9523 0.8241 0.8992 0.8856 

/+(1.16) 

1.0047(123) 

1.0000 0.7892 0.8900 0.9530 

/o(l) 

0.9026(072) 

1.0000 0.9650 0.8682 

/o(1.08) 

0.8609(077) 

1.0000 0.9519 

/o(1.16) 

0.8254(094) 

1.0000 


therefore negligible. Although this method varies the light valence- and sea-quark masses 
together, the shifts are primarily due to the different valence-quark mass. 


V. DETERMINATION OF \Vcb\ 

A. Synthetic data 

The preferred chiral-continuum £t results for h^{w) and h-{w) are continuous functions 
of w at zero lattice spacing and physical quark masses. Via Eqs. (2.5) and (2.6), we can 
obtain the corresponding functions for f+{w) and fo{w). As discussed above, the errors are 
under control for w < 1.2, i.e., where we have lattice measurements. Following Refs. [14, 25], 
we proceed to extend our results to the full kinematic range by generating synthetic data 
f+{wj) and foiwj) for a hnite set of w values, Wj. Because the functions are described by 
only six independent functions (in the physical limit), we can only generate six such data 
points. Generating more would just lead to a covariance matrix of low rank. We choose the 
values of Wj = 1, 1.08, and 1.16, for /+ and /o, to cover the kinematic range of the lattice- 
QCD calculation. The values of /+(wj) and fo{wj), as well as the matrix of correlations 
among them, are given in Table VII. 


B. 2; expansion 

Experimental measurements of the form factor are available over a larger kinematic range 
of w [1, 1.58] than the lattice values [1,1.16], but experimental errors are largest where lattice 
errors are small and vice versa. Although the value of /+ at a single tc-value suffices for 
obtaining iWbl, a better strategy is to fit both sets of data simultaneously to a common 
htting function in which \ Vcb\ is a free parameter that multiplies all the lattice values and is 
determined in the fit [14, 25]. This approach minimizes the uncertainty in \ Vcb\ by combining 
all of the available experimental and lattice information. Further, a comparison of the shapes 
of the experimental and lattice results as a function of w provides a valuable consistency 
check that is not available when using only a single recoil point. 

For this purpose we need a model-independent parameterization to carry out the nec¬ 
essary interpolation/extrapolation. The 2 ; expansion of Boyd, Grinstein and Lebed (BGL) 
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[11] is just such a parameterization. It builds in constraints from analyticity and unitarity. 
It is based on the conformal map 


z(w) 


a/ 1 + w — a/2 
a/1 + w + y/2 ’ 


(5.1) 


which takes the physical region w G [1,1.59] to ^ G [0,0.0644]. It pushes poles and branch 
cuts relatively far away to \z\ ~ 1. Form factors are then parameterized as 


fiiz) 


1 

Pi{z)(j)i{z) 


OO 



n=0 


(5.2) 


where the Pi{z) are the “Blaschke factors” containing explicit poles {e.g., a or B* meson) 
in the channel variable and the 0* are the “outer functions”, whose purpose is described 
below. The only unknown parameters are the polynomial coefficients ai^n- In this work, we 
do not introduce any pole, so Pi{z) = 1.^ The choice of outer functions is arbitrary as long 
as they are analytic functions that do not introduce poles or branch cuts; the 0i just affect 
the numerical values of the series coefficients, Uj. For /+ and /o, we use 

(l)+{z) = <I)+(1 + zf{l - zY^‘^[{l + r)(l - z) + 2y/r{l + z)]~^ , (5.3) 

(j)o{z) = <I)o(l + ^)(1 - zf/\l + r)(l - z) + 2y/r{l + z)]~^ , (5.4) 

such that, numerically, $0 = 0.5299 and <I>+ = 1.1213 [11]. With this choice, the bound on 
the series coefficients from unitarity takes a particularly simple form: 


N 

\ai,n? < 1 , 

n=0 


(5.5) 


where this bound holds for any N . This bound, in combination with the small range of \z\, 
ensures that only a small number of coefficients is needed to parameterize the form factors 
over the entire kinematic range to high precision. 

To implement the z expansion, we start from the synthetic data for /_|_ and /o at 2 ; values 
corresponding to wj = 1, 1.08, and 1.16, choose a truncation N and £t to determine the 
coefficients for n = 0,..., N. These coefficients are then used to parameterize the form 
factors over the full kinematic range. We hnd we need only the hrst few coefficients in the 
expansion to obtain a stable fit with a good p value. The kinematic constraint requires 
/+ = /o at = 0 where 2 ; ~ 0.0644. It is interesting to £t the data without the constraint 
to see to what extent it is automatically satished. The result for = 3 in the left panel 
of Fig. 8 shows that the data satisfy the constraint much better than our statistics would 
suggest. Nonetheless, in subsequent hts, we include the constraint to reduce the form-factor 
errors at g^ = 0. The constraint is imposed by expressing the parameter ao,o in Eq. (5.2) in 
terms of the other series coefficients. Table VIII shows the series coefficients and goodness- 
of-£t obtained for hts of the lattice form-factor data imposing the kinematic constraint with 
N = 2-4. For the hts at cubic and quartic order in the z expansion, we have more parameters 
than data, but the unitarity bound in Eq. (5.5) justihes imposing a prior with central value 
0 and width 1 on the coefficient(s) of the cubic (and quartic) term(s). 


^ We have checked that including a pole located at the theoretically-predicted B* mass [43] does not 
appreciably change the z-&t result. 
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FIG. 8. Result of the z-expansion fit of the lattice form-factor values without (left) and with (right) 
the kinematic constraint f+{q‘^ = 0) = foiq^ = 0). The expansion is truncated after the cubic term. 
The solid error band is for /+, while the slashed band is for /q. Without imposing the constraint, 
we find that it is nonetheless satisfied to a high accuracy. 


TABLE VIII. Coefficients of the z expansion for hts to the lattice form factors including the 
kinematic constraint f+{q‘^ = 0) = fo{q‘^ = 0). For completeness, the inferred value and error in 
ao,o is quoted. We also show the zero-recoil form factor ^(1). The results for different truncations 
N are virtually identical. The unusnally low (augmented) comes about because these fits 
essentially behave like solves. This happens because the kinematic constraint is so nearly perfectly 
satisfied already at the quadratic level, N = 2. Higher-order terms with N = 3 and 4 provide no 
further improvement and, hence, no change. 



N = 2 

N = 3 

iV = 4 

O'+,0 

0.01262(10) 

0.01262(10) 

0.01262(10) 

O'+A 

-0.097(3) 

-0.097(3) 

-0.097(3) 


0.50(14) 

0.50(17) 

0.50(17) 


- 

-0.06(90) 

-0.06(90) 


- 

- 

-O.O(l.O) 

ao,o 

0.01142(14) 

0.01142(14) 

0.01142(10) 

00,1 

-0.060(3) 

-0.060(3) 

-0.060(3) 

00,2 

0.31(15) 

0.31(15) 

0.31(15) 

Oo,3 

- 

0.06(91) 

0.06(91) 

00,4 

- 

- 

O.O(l.O) 

^(1) 

1.0541(83) 

1.0541(83) 

1.0541(83) 


0.1/1 

0.0/1 

0.0/1 
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TABLE IX. Central values, errors, and correlation matrix for the parameters of the cubic fit to /+ 
and /o including the kinematic constraint at = 0. 



value 

Correlation matrix 

“- l -,0 “+,1 “+,2 “+,3 “ 0,1 “ 0,2 “ 0,3 


0.01262(10) 

1.00000 0.21726 0.07203 0.00387 0.19347 0.15590 -0.00364 


-0.0969(34) 

1.00000 -0.47505 0.25544 0.80946 -0.26302 -0.18212 


0.50(17) 

1.00000 -0.45415 -0.43845 0.85491 0.25116 


-0.06(90) 

1.00000 0.11415 -0.15582 0.21768 

“ 0,1 

-0.0597(29) 

1.00000 -0.42932 -0.03556 

“ 0,2 

0.31(15) 

1.00000 -0.06062 

“ 0,3 

0.06(91) 

1.00000 


The truncation of the ^ expansion introduces a possible systematic error. We take this into 
account by increasing the truncation order until the central values and errors stabilize. At 
this point, the errors from the fit reflect the truncation error, and do not need to be counted 
separately. Table VIII shows that the £t has stabilized by quadratic order. We therefore 
take the cubic fit, shown in the right panel of Fig. 8, as our preferred parameterization. 
Table IX gives the central values, errors, and normalized correlation matrix for the series 
coefficients a*. This information can be used to reproduce our results for f+{w) and fo{w) 
over the full kinematic range, and in particular, in combined lattice-and-experiment fits to 
obtain |Wb|- 

We compare our form-factor results with those of the most recent lattice-QCD calculation 
of 5 —?• Diu at nonzero recoil in Fig. 9. Although this earlier calculation was performed 
in quenched QCD, and thus is subject to an unquantifiable systematic due to the omission 
of sea-quark effects, it uses step-scaling [44] to control heavy-quark discretization effects, 
plus multiple light-quark masses and lattice spacings to control the mild chiral-continuum 
extrapolation [12]. Thus it is the best calculation so far for B —)■ Diu at nonzero recoil. The 
two calculations agree for all w values, although the slope of f+{z) is somewhat steeper for 
the (2-1-1)-flavor result reported here. 


C. Determination of \Vcb\ 

To obtain iWbl, we need lattice results for the form factors and experimental values for 
VEw\Vcb\f+{w). Because the experimental value of the form factor at zero recoil suffers 
from kinematic suppression, we prefer to fit the theoretical and experimental data over the 
entire kinematic range. For this work, we use the 2009 i?-tagged data from the BaBar 
collaboration [8], because it is the most precise to date.^ Reference [8] reports a systematic 
error of 3.3% at small w. For present purposes, we take 3.3% over the entire kinematic 
range with 100% correlation and combine this systematic in quadrature with the reported 
(uncorrelated) statistical errors [45]. 

Although the BaBar collaboration has applied some radiative corrections to their pub¬ 
lished data, additional electroweak effects still remain. These include a Sirlin factor for the 

^ The Belle experiment presented preliminary measurements of l?_Ew|Fcf,|/+(w) at ICHEP 2014 [9]. Once 
these are finalized, our form-factor coefficients from Table IX can be used to update iRbl from a joint 
lattice-experiment fit with both the Belle and BaBar data (including experimental correlations). 
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FIG. 9. Comparison of lattice-QCD results for the B —)• Diiy form factor G{z) at nonzero recoil 
from this work (curves with error bands) and Ref. [12] (points with error bars). Errors on the data 
points from Ref. [12] include all uncertainties except for the unquantifiable error due to omitting 
sea-quark effects. 


VFq and WZ box diagrams [46] and a further Coulomb correction for final-state interactions 
in decays. The BaBar collaboration reports that 37% of the decays in their data sample 
were B^s, which results in a QED correction factor in the amplitude of 1 -|- O.STa/(27r). We 
have assigned an uncertainty of ±0.005 to this correction to account for omitted electromag¬ 
netic effects at intermediate distances. When combined with the Sirlin factor rjEw = 1-00662 
the net electroweak correction becomes ?)ew = 1-011(5). (We prefer to use Q{w) to denote 
the purely hadronic form factor, so in our notation f]Ew\Vcb\Q{w) corresponds to the quan¬ 
tity often reported as \ Vcb\G{w), and the ratio of experimental to theoretical values must be 
divided by fjEw to get |Wb|-) 

Before performing a joint £t to the lattice and experimental data, we compare the values 
of the shape parameters to check for consistency. The left panel of Fig. 10 plots the l-cr 
constraints on the curvature a+, 2 /a+,o versus slope a+q/a+,o obtained from separate N = 3 
z-expansion fits of the lattice data and the 2009 BaBar experimental data. The results are 
consistent, but the lattice data constrains the shape much better: this is both because the 
lattice points are very precise at low recoil, and because they are more correlated between 
w values. Given this consistency, we now proceed with the determination of \Vcb\ from a 
combined fit of the two data sets. 

Table X shows the series coefficients and goodness-of-£t obtained for combined hts of the 
lattice and experimental data, imposing the kinematic constraint, for N = 2-4. Again, the 
£t, and in particular the error on iWbl, stabilizes by quadratic order. We choose X = 3 for 
our preferred £t, and plot the result in Fig. 11. 
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FIG. 10. Left: One sigma contour plots showing the correlation between the normalized slope 
a+^i/a+^o and normalized curvature a+^ 2 / 0+,0 from N = 3 z-expansion fits to either the BaBar 
experimental data alone, our lattice QCD results alone, and a joint fit to both. Right: vector form 
factor /+ obtained from separate z-expansion fits of the 2009 BaBar experimental data (hatched 
band) and lattice form factors (solid band). 


TABLE X. Best-fit values of the z-expansion parameters for different truncations N from a joint 
fit to experimental data and lattice values. For completeness, the inferred value and error in ao,o 
is quoted. We also show the zero-recoil form factor ^(1) and |I46|- 



N = 2 

N = 3 

iV = 4 

o-+,o 

0.01260(10) 

0.01261(10) 

0.01261(10) 

O'+A 

-0.096(3) 

-0.096(3) 

-0.096(3) 

0-+,2 

0.37(8) 

0.37(11) 

0.37(11) 


- 

-0.05(90) 

-0.05(90) 


- 

- 

-O.O(l.O) 

00,0 

0.01140(9) 

0.01140(9) 

0.01140(9) 

Oo,l 

-0.059(3) 

-0.059(3) 

-0.059(3) 

O0,2 

0.18(9) 

0.19(10) 

0.19(10) 

00,3 

- 

-0.3(9) 

-0.3(9) 

00,4 

- 

- 

-O.O(l.O) 

^(1) 

1.0527(82) 

1.0528(82) 

1.0528(82) 


0.0396(17) 

0.0396(17) 

0.0396(17) 

X^/df 

8.4/10 

8.3/10 

8.3/10 
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FIG. 11. Result of the preferred joint fit of the BaBar experimental data together with the lattice 
form factors. The plotted experimental points have been divided by our best-fit value of f/Ew|f4fe| 
and converted to /+. 


D. Comment on the CLN parameterization 

The standard approach used by experimentalists to obtain \ Vcb\ is to use the Caprini, Lel- 
louch, Neubert (CLN) parameterization [10] to extrapolate the experimental data to w = 1. 
Caprini, Lellouch, and Neubert use heavy-quark symmetry to derive more stringent con¬ 
straints on the coefficients of the ^-parameterization through 0{z^), resulting in a function 
with only two free parameters, /+(0) and pf: 

= 1 - 8plz + {blpl - 10)^2 - (252p2 - 84)^=". (5.6) 

/+(b) 

Use of the CLN parameterization in our analysis does not reduce the quoted errors in \Vcb\ 
despite the introduction of additional theoretical information. 

The numerical values of the coefficients in Eq. (5.6) have theoretical uncertainties which 
can be estimated from the information given in tables and plots from Ref. [10]. To the best 
of our knowledge, however, CLN fits to experimental data do not incorporate the theoretical 
uncertainties discussed in Ref. [10], and may therefore be underestimating the uncertainty 
in |Ucfe|- We have attempted to quantify the uncertainty from the use of the CLN form by 
incorporating the theoretical uncertainties in the CLN parameters via Bayesian priors. We 
did not find any difference in the error on \ Vcb\ obtained from fits with and without including 
these theoretical uncertainties at the current level of precision. This is primarily because 
the B —)■ Diu data displays little evidence of curvature in 2 ; within the present errors, and 
does not constrain the coefficient of the z^ term. Nevertheless, we do not quote the results 
of our CLN fits in this work because we are more confident in the errors obtained from the 
model-independent ^^-parameterization, Eq. (5.2), which can be used to obtain \Vcb\ even as 
the experimental and lattice uncertainties become arbitrarily more precise. 


29 











VI. DISCUSSION AND OUTLOOK 


We obtain 

iWbl = (39.6 ±1. 

^QCD+exp ± 0.2qed) X 10 ^ (6.1) 

from our analysis of the exclusive decay B Dlv aX nonzero recoil, where the hrst error 
combines systematic and statistical errors from both experiment and theory and the second 
comes from the uncertainty in the correction for the final state Coulomb interaction in the 
decays. Because we provide the series coefficients of a z parameterization and their 
correlations, the result for \Vcb\ in Eq. (6.1) can be updated whenever new experimental 
information becomes available. 

The combined error from lattice and experiment in \Vcb\ is about 4%. Because this error 
is obtained from a joint z-fit, the theory and experimental errors cannot be strictly disen¬ 
tangled, but they can be estimated as follows. In the right panel of Fig. 10 we plot the 
determinations of /+ from separate z fits to the lattice form factors and to the experimental 
data. Inspection of the error bands shows that the combined error, which determines the 
uncertainty on iWbl, is smallest at about z 0.025 {w ~ 1.2). At this point, the experi¬ 
mental error is about 3.9% and the lattice error is about 1.4%. (Note that combining them 
in quadrature yields a total that is close to the 4% lattice-f-experiment error on \Vcb\ from 
the joint fit.) Thus the experimental error currently limits the precision on \Vcb\ from this 
approach. The dominant uncertainty in the experimental data is the assumed 3.3% system¬ 
atic error, which is used for all w values in the joint £t. Now that lattice-QCD results for 
the B —)■ D^v form factors are available at nonzero recoil, however, it is clearly worthwhile 
to study and improve the systematic errors in the experimental data at medium and large 
recoil. 

It is interesting to compare the above nonzero-recoil result with the result based on the 
standard method that uses only the zero-recoil extrapolation of the experimental and theo¬ 
retical form factors. The 2 ; expansion fit to lattice-only data gives ^(1) = 1.054(4)stat(8)syst- 
The BaBar collaboration quotes f]Ew\ycb\Q{X) = 0.0430(19)stat(14)syst [8] from its i?-tagged 
data, which gives \Vcb\ = (40.8 ± 0.3qcd ± 2.2exp ± 0.2qed) x 10“^. The result is con¬ 
sistent with the value from nonzero recoil, but the error is larger, as expected. Our 
zero-recoil form factor is consistent with a previous, preliminary Fermilab/MILC result 
of ^(1) = 1.074(18)stat(16)syst [13], but with significantly smaller uncertainties due to the 
use of a much larger data set with several lattice spacings and lighter pious. We also note 
that the systematic error estimate for the earlier result did not include an estimate of the 
heavy-quark discretization errors, one of the larger contributions to the error in our new 
result. 

We compare our result for \Vcb\ with other published determinations from inclusive and 
exclusive decays in Fig. 12. Our result is consistent with the determination from our compan¬ 
ion analysis of 5 —?■ D*ii' at zero recoil, \Vcb\ = (39.04 ±0.53qcd ±0.49exp ± 0.19 qed) x 10“^ 
[5]. The errors on \ Vcb\ from the current work are larger, however, because of the larger errors 
in the experimental data. Our result is 1.5cr lower than a recent inclusive (non-lattice) de¬ 
termination, \Vcb\ = (42.4±0.9thy+exp) X 10“^ [6], which is also based on several experiments 
and employs data at nonzero recoil. 

We also plot the result for \Vcb\ in Fig. 12 determined from only our zero-recoil lattice 
data, but using the best experimental knowledge of the extrapolated quantity f]Ew\ycb\Q{X)- 
The HFAG average value W;,|^(l) is 0.04264(72)stat(135)syst [3], which combines five 
experimental measurements from ALEPH [47], Belle [48], BaBar [8, 49], and CLEO [50]. 
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Fermilab/MILC ’ 15 + HFAG ’\A,B^D,w=\ 

Fermilab/MILC ’ 14 + HFAG 'U,B^D,w=\ 

Gambino & Schwanda ’ 13, B —> X inclusive 
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FIG. 12. Comparison of exclusive and inclusive determinations of \ Vcb\ x 10^. Triangles denote an 
extrapolation to zero recoil, while squares use data over a wide kinematic range. The color code 
is black, blue (dark gray), and orange (light gray) for B —)• Diu, B —>■ D*iv, and B —)■ Xciv, 
respectively. 
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From this value we obtain \Vcb\ = (40.0 ± O.Sqcd ± l.dgxp ± 0.2 qed) x 10“^. This error is 
smaller than that from the analysis at nonzero recoil, thanks to the additional experimental 
information, but only by about 10%. Thus combining lattice data at nonzero recoil with a 
single experiment reduces the error on \Vcb\ by almost as much as adding zero-recoil data 
from several experiments. Clearly the error on \ Vcb\ from B -A- Diu at nonzero recoil can be 
further reduced via a joint fit of the lattice form-factor data with additional experimental 
measurements once correlations are available. 

An interesting byproduct of our combined ^-expansion £t to obtain \Vcb\ is an improved 
determination of the B ^ D form factors f+iq"^) and foiq"^)- Because the lattice form factors 
are most accurate at high q‘^, while the experimental measurements are most accurate at 
low they provide complimentary constraints on the form-factor shape. Table XI provides 
the 2 ;-fit coefficients and correlation matrix from our preferred combined lattice-experiment 
fit used to obtain our result for \Vcb\ quoted in Eq. (6.1). These represent our current best 
knowledge of f+iq"^) and foiq"^) ior B ^ D semileptonic decays, and can be used in other 
phenomenological applications. Here we use the results in Table XI to update our calculation 
of the ratio B{B ^ Dti>)/B{B —)■ Diu) in the Standard Model [15]. We obtain 

i?(Z9) = 0.299(11), (6.2) 

which agrees with our previous determination R{D) = 0.316(12)(7) in [15], but is 2.0a 
lower than the BaBar measurement R{D) = 0.440(58)(42) [51]. The error in our new 
determination of R{D) is about 20% smaller than in Ref. [15], primarily due to the inclusion 
of the experimental information on the shape of /+ from the joint z-fit. 

The dominant errors in the lattice form factors come from statistics, matching, and the 
chiral-continuum extrapolation, and can be reduced through simulations at smaller lattice 
spacings and at physical quark masses and from further study of the matching factors. The 
MILC Collaboration is currently generating (2-|-1-1-1)-flavor HISQ ensembles with physical 
light quarks [52], which we anticipate using for future calculations of i? —)■ D^*'^ form fac¬ 
tors. Heavy-quark discretization errors are also important. They can be reduced with a 
more improved heavy-quark action such as that proposed in Ref. [53], and work on this is 
underway [54, 55]. 


31 







TABLE XL Central values, errors, and correlation matrix for the parameters of the joint cubic fit 
to the synthetic lattice data for /+ and /o (including the kinematic constraint at (f' = 0) and the 
experimental measurements of fiEw\ycb\f+{w)- 



value 

Correlation matrix 

O +,0 0+4 a +4 0+4 00,1 Oo ,2 Oo,3 


0.01261(10) 

1.00000 0.24419 -0.08658 0.01207 0.23370 0.03838 -0.05639 


-0.0963(33) 

1.00000 -0.57339 0.25749 0.80558 -0.25493 -0.15014 


0.37(11) 

1.00000 -0.64492 -0.44966 0.66213 0.05120 


-0.05(90) 

1.00000 0.11311 -0.20100 0.23714 

00,1 

-0.0590(28) 

1.00000 -0.44352 0.02485 

00,2 

0.19(10) 

1.00000 -0.46248 

00,3 

-0.03(87) 

1.00000 
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Appendix A: Heavy-quark mass correction 


Heavy-quark masses {k values) are determined by requiring that the kinetic masses of 
the Dg and Bg match their experimental values. The three-point and two-point functions 
in this study were computed with k, values from a preliminary tuning. Final tuned values 
differed slightly [5], as shown in Tables II and V. 

We therefore need to adjust the form factors and w values accordingly. This is done 
by repeating the computation of the ratios R+, Q+{p), and R^{p) on the a ~ 0.12 fm, 
m' = Q.2m'g ensemble for a few values of and in the vicinity of the desired, tuned 
values. These results permit calculating the derivatives of the form factors with respect to 
the quark masses. We assume that these results, expressed in dimensionless terms, can then 
be used to adjust form factors in our other ensembles. 

From Eqs. (2.18)-(2.20), we see that we have the option of computing and applying 
these adjustments before or after matching with the p factors. Because of the simplifying 
steps taken in Appendix B for pvi{w)/pv^iX) and pvi{w)/pv^iw), we choose to make the 
adjustments directly on unmatched quantities. From Eqs. (2.22) and (2.23), one sees that 
it is convenient to study the mass dependence of 




a/ i?+Q+, 
R- Xf 


(Al) 

(A2) 


and x'jS-. 

Heavy-quark symmetry suggests that interpolations in inverse quark masses will imple¬ 
ment the quark-mass tuning most smoothly. With the Fermilab method [16], the quark 
mass is identihed with the kinetic mass: 


1 _ 2 1 
77120 7?7oa(2 + TTloO) TTloO + 1 ’ 


(A3) 


where we compute the bare quark mass moa from the tadpole-improved, tree-level formula 


TUqU = 


Uq \2k 



(A4) 


Here, uq is the tadpole parameter, and is the value of k such that the lightest pseudoscalar 
meson mass vanishes. Thus, below we compute slopes of S'+, S'-, and aJjS'_ with respect to 
fy = 1/(7T72 c’"i) and fy = l/(7T72fe7"i). The results of the computations with varying quark 
masses are shown in Figs. 13 and 14. 

Because the corrections in the charm and bottom masses are small, it suffices to work to 
hrst order in the inverse mass shift. Heavy-quark symmetry also suggests that the leading 
mass dependence of S'+(7n = 1) = y/R+ is quadratic, of the form (fy — ^b)^. Therefore, 
the leading shift in fy is suppressed by fy, and the leading shift in fy is suppressed by fy. 
Below, we neglect the former effect but keep the latter, since it is suppressed only by fy. 
Furthermore, by construction x'^jS- —)■ 0 as tc —)■ 1 for all quark masses, and, therefore, the 
derivative with respect to fy also vanishes at w = 1. On the other hand, neither S- nor its 
derivatives vanish aX w = 1. 
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FIG. 13. Heavy quark mass dependence on the a ~ 0.12 fm, rh' = 0.2mg ensemble at momentum 
27r(l, 1,0)/L. Left to right: 5+, S-, and x'jS-, respectively, vs. inverse charm-quark kinetic 
mass 
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FIG. 14. Heavy quark mass dependence on the a ~ 0.12 fm, rh' = 0.2m^ ensemble at momentum 
27r(l, 1,0)/L. S+ (left) and S- (right) vs. inverse bottom-quark kinetic mass = {m 2 bri)~^. 


Because of the narrow range of w, 1 < w < 1.16, for our data, one should expect a linear 
approximation in w to suffice for the quark-mass adjustments. Indeed, only requires 

a quadratic, as shown in Figs. 15 and 16. Therefore, we introduce 


dS+ 

dS_ 

d{x)S.) 

die 

dw 

Wc 

dS+ 

dib 

dS_ 

dib 


r+,iAw - 1) , 

(A5) 

r_,o,c + r-xc{w - 1), 

(A6) 

rx,iA'^ - 1) + ’^x,2,c(w - 1)^ , 

(A7) 

rw,iA^ - 1) 1 

(AS) 

r+,o,b + r+^iAi^ - 1) 1 

(A9) 

r_,o,b + r-^iAw - 1). 

(AlO) 


(The notation for the slope parameters rf^n,q encodes a form factor label /, a polynomial 
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FIG. 15. Effect of heavy-quark mass shifts on the a ~ 0.12 fm, m' = 0.2m^ ensemble as the 
charm-quark mass parameter is increased from Kc = 0.1254 to 0.1280. Left to right: dS^, dS-, 
and d{x‘jS-), respectively, vs. w — 1. 




FIG. 16. Effect of heavy-quark mass shifts on the a ~ 0.12 fm, m' = 0.2m^ ensemble as the 
bottom-quark mass parameter is increased from = 0.0860 to 0.0901 vs. w — 1. Left to right: 
dS+ and dS-, respectively, vs. w — 1. 


coefficient index n, and a quark mass label q.) Fits to our data then yield 


r+,1,, = -0.72(5), (All) 

r_,o,c = 0.102(11), (A12) 

= -0.23(14), (A13) 

= 0.851(14), (A14) 

r,,2,e = -1.22(16), (A15) 

r+,o.f, = 0.0042(41), (A16) 

= 0.07(21), (A17) 

= -0.49(5), (A18) 

= -0.46(69), (A19) 


As discussed above, we expect to be of order ^c-, or approximately 0.83. In fact, it is 
much smaller. 

We compute the correlation functions at discrete values of the recoil momentum of the D 
meson, resulting in discrete values Wi, which are determined from Eqs. (2.15) and (2.17). The 
recoil variable Wi is determined dynamically from diagonal vector current matrix elements 
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FIG. 17. Effect of heavy-quark mass shifts on the a ~ 0.12 fm, m' = 0.2m^ ensemble. Shift in 
recoil variable Wi vs. w — 1 as the charm-quark mass parameter is increased from Kc = 0.1254 
to 0.1280. 


TABLE XII. Unadjusted and adjusted values of w, h+, and h- for the a ^ 0.12 fm, m' = 0.14m'^ 
ensemble. Eor this illustration only statistical errors are shown. 


momentum 

raw w 

shifted w 

raw /i+ 

tuned h+ 

raw h- 

tuned h- 

000 

1 

1 

1.0391(53) 

1.0390(53) 

- 

- 

100 

1.0465 

1.0426 

0.9812(61) 

0.9849(62) 

0.0041(111) 

0.0008(111) 

110 

1.0896 

1.0822 

0.9388(70) 

0.9457(70) 

0.0072(134) 

0.0031(135) 

111 

1.1299 

1.1191 

0.8978(97) 

0.9073(97) 

0.0139(165) 

0.0092(166) 

200 

1.1553 

1.1424 

0.8789(120) 

0.8898(121) 

0.0336(205) 

0.0286(206) 


involving the D meson, so it varies with the charm quark mass, but not the bottom quark 
mass. We take the convention that when we shift both quark masses, we shift Wi to w[ and 
we shift Sj^{wi) to and similarly for S- and x'jS-. As can be seen from Fig. 17, the 

data support a linear approximation for the shift in Wi also. 

The effect of the kappa adjustment on the values of w, and is illustrated in 

Table XII for the a ~ 0.12 fm, rh' = 0.14m'^ ensemble where the tuning adjustment decreases 
both Kc and ki, from their simulation values. 
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Appendix B: Heavy-quark discretization effects 


We use the heavy-quark effective theory (HQET) to derive the form of heavy-quark dis¬ 
cretization effects [18, 39]. In this appendix, we apply this formalism to derive the matching 
procedure from lattice gauge theory to continuum QCD, cf. Eqs. (2.18)-(2.20). We also use 
it to derive power-law discretization effects, both at nonzero recoil {w > 1) and at zero recoil 
(tc = 1) where heavy-quark symmetry is more constraining. In the last subsection of the 
appendix, we also present numerical estimates for the discretization errors. 


1. Formalism 


We describe the underlying lattice gauge theory (LGT) with an effective Lagrangian, 
asserting the relation 


. - hD\h hs ■ Bh h[D°[,iEa\h hsaii{DJ^,iE^}h 

Clgt = h{iv ■ D - mi)h + ^ ^ ^ + 


2m9 


2m 


B 


Smjj 




(Bl) 


where = can be read “has the same matrix elements as.” Here, n is a four vector specifying 
the rest frame of the heavy-light meson, such that = —1; the heavy-quark held h satishes 
'fih = ih] and = —iaapl2. Then, = D^ + v^ v-D is the covariant derivative orthogonal 

to n, = (5" + is the chromomagnetic held (in the v frame), and 

Ef^ = -VaE^^ is the chromoelectric held (in the v frame). The HQET description for 
continuum QCD has the same structure 


, , h zuhs 

m)h+ ^ + 


Bh 


2m 


2m 


^ ZDh[D1,iEa]h ^ ZEhSais{D1,iE^}h ^ 




Am"^ 


— h{iv ■ D 

/'in /'III /\'iri^ 0.111“ 

(B2) 

In this framework, matching and improvement boil down to choosing the parameters of the 
lattice Lagrangian, such that the Eq. (Bl) reproduces Eq. (B2) term by term. 

The rest mass mi does not inhuence matrix elements or mass splittings. In the Fermilab 
method, therefore, one tunes k so that 


1 _ 1 
2 m 2 2m ’ 


(B3) 


and csw so that 

— = ^ = (B4) 

2m B 2m 2m 

where the second equality follows because = 1 + 0{as)- In this work, we tune k via 

the heavy-strange meson mass; for details of our procedures, see Appendix C of Ref. [5]. 
Furthermore, we choose csw at the tadpole-improved tree level, which makes the coefficient 
of the 0{as) error in Eq. (B4) small [56]. 

The Fermilab vector current, Eq. (2.8), has an HQET description too. Through dimension 
four [18] 


14 


V'‘ = Cy°'^v'‘c„:K + ^ h — 


(B5) 


a=l 
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The continuum-QCD current can be described in the same way albeit with different C 
and B coefficients, denoted in this paper by omitting the label “LGT.” Then 
if the Z factors are chosen to be [18] 

Zy^^=CyJC]f^^, (B6) 

Zy. = Zy^=CyJC]ff^. (B7) 

In practice, of course, such matching is only approximate. For example, the one-loop calcu¬ 
lation of py4, explained in Sec. HIE, leads to a matching error of order 

With the Fermilab currents most of the fourteen dimension-four coefficients vanish 

at the tree level; the same holds for continuum QCD and the By a- The exceptions multiply 
the operators 


Qyi = 

= Cyi'y^lpxK, 
Qv 4 = -v'^Cylpx'by, 

At the tree level, their coefficients are 


7 rLGT _ y rLGT _ 
— ^V±-Dv2 = 


7 rLGT _ y rLGT _ 
^V«^y4 — ^Vj_-Dy^ = 


1 


2m3b‘ 

1 


2m^c 

The improvement terms in the current, namely di, are chosen so that 

+ O^agCi), 


(B8) 

(B9) 

(BIO) 

(Bll) 

(B12) 

(B13) 

(B14) 


2m3 2m2 

for operators with label a G {1, 2,4, 5}. The other By^^'^ are of order Og from the outset. 


2. Matching factors 

Equation (2.18) is well known from earlier work [18, 19]. To establish Eqs. (2.19) 
and (2.20), let us start by dehning h^^{w) for the lattice current in analogy with 
the decomposition in Eq. (2.4). These form factors h^^{w) are not the right-hand sides of 
Eqs. (2.22) and (2.23). The task here is to show how the ratios cancel some of the cutoff 
effects in /^[^^"'"(w). Sometimes it is convenient to choose arbitrary v and n' when working 
out consequences of the HQET. The kinematics of our lattice-QCD correlators correspond 
to n = {i, 0) and v' = {iw, v'). 

The simplest case is the dehnition of the velocity via D{0) —> D{p) matrix elements: 

" (to + l)ftyT(p) _ („ _ l)gLGT(p( ,„_il. 1““^ 

because = 0 for a flavor-conserving transition. This property follows from time-reversal 
invariance of the chosen current and arises independent of any matching considerations. The 
expression for w in Eq. (2.17) then follows immediately from w'^ = 1 + (when n = 0). 
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Similarly, the other ratios are 


Q+{p) = 


{w + — (w — l)h^^'^{w) 


R-ip) = v' 


^LGT(y^) _ 


{w + — (w — l)h^^'^(w) ’ 


(B16) 

(B17) 


with w = w(p). These form factors, of course, are for the flavor-changing process. 

Using the trace formalism explained in Ref. [39], it is straightforward to obtain the 
following expressions for h^^'^(w): 


h^^^(w) = C^^^(w)S(w) + [2e3(wi) - MH] - , 

(B18) 

h^y^iw) = liw + l)C^^^iw)E{w) + [2e3(w) - Ae(rc)] - B'!:^^{w)M{w), 

(B19) 


neglecting higher-dimension terms. The leading-dimension, short-distance coefficients are 

+ liw - l)C^^yw), (B20) 

Cy^{w) = - Clf^^iw) -{w + (B21) 


The contain several of the fourteen By^ in Eq. (B5), and the detailed ex¬ 

pressions are not illuminating. The Isgur-Wise function yw) and its generalizations ^ 3 ( 10 ) 
and 

E{w) = yw) + E 2 Ai{w) + Eb [3v43(m;) -h 2 (m; - l)R 2 (tf)] (B22) 

parameterize the long-distance physics. In the context of lattice gauge theory, their dis¬ 
cretization effects arise only from the light degrees of freedom. In Eq. (B22), ,^(1) = 1 and 
Ri(l) = ^ 3 ( 1 ) = 0 by flavor conservation in HQET. In order to have compact formulas, the 
function H contains some short-distance information, namely the mass combinations 


Ux = 


1 

2 mxc 


1 

2 mxb' 


X e { 2 , B, 3 }, 


(B23) 


which depend on the short-distances a and nigh 

When using HQET to describe the heavy-quark limit of continuum QCD, the algebra is 
identical. The difference lies in the short-distance coefficients: in the notation used here, 
Cy^'^ etc. simply lose the superscript “LGT”. Further, discretization effects of the light 
degrees of freedom disappear from the HQET quantities A, ^(w), ^ 3 ( 10 ), and Ai{w). 

To derive the matching factors, we focus on the leading-dimension term. Then one hnds 


Q+{p) 


R~{p) 


v' Cyf^{w) -1- (w -1- 

tc -f-1 C\y'^{w) 


(B24) 

(B25) 
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Thus, to match these quantities to continuum QCD, one must multiply Q+ and R- by 


Pv*{w) _ Pv^iiw) _ Cv^^{w) 

Pv<l) ^ pv,il) ~ C^^^iw) Cv,il) ’ 

Pv^jw) ^ P\/,(w) _ Cy^jw) + (w + l)(gt,/(w) 

Pv4w) ~ pV||(w) ^ Cv^^{w) 


(B26) 

(B27) 


respectively, to obtain Q+ and IZ^ in Eqs. (2.19) and (2.20). 

One-loop calculations of the w dependence of these coefficients are not available, however. 
(The algebra with p 7 ^ 0 is much more voluminous.) We shall proceed with a further 
approximation for each of the two factors multiplying Q+ and i?_. By construction in 
Eq. (B26), 

PVn (U!) 

^^ = l + Oiasiw-l)). (B28) 

PV||(1) 

Because the w dependence arises only from the vertex diagram—the tadpoles on the legs 
cancel—the computed coefficient should, like those in Table V, be small. Furthermore we 
note that w — 1 < 0.16 and that the w dependence disappears when rrica —)■ 0. Hence we 
neglect this one-loop contribution and take py^^w)/pv4(l) = 1. For the 0{as) error we use 
the following form: 

= 1 ± ay{2/- l)m 2 ca, (B29) 

where the values for ay {2/a) are listed in Table V, and 


= (B30) 

is an upper bound on the size of the observed one-loop corrections to pyA{l). In the mass 
region of interest, < Pyl^ax- 

Equation (B29) gives an estimate of the error in the ratio py4{w)/pyA{l). The zero- 
recoil py 4 (l) is calculated at one-loop order in lattice perturbation theory and tabulated in 
Table V. We estimate the 0(a‘^) truncation error, in the spirit of Ref. [5], by taking the 
coefficient as twice the largest hrst-order coefficient, = 0.2. Hence, the error due to 

omitted higher order corrections is estimated as 

± 2pb„„“?'(2/o). (B31) 

The two errors are combined in quadrature to obtain the total systematic error in pyA[w): 


± pyA{l)^[p^y\^^^ay{2/a){w - l)m2ca]2 + [ 2 pW ^^^a^(2/a)/py4(l)]2 . (B32) 

For the factor in Eq. (B27) for i?_, note that most of our ensembles have rrica < 0.4 and 
recall that as rrica —)■ 0 with mi,a hxed, the short-distance coefficients of the HQET with two 
heavy-quark helds tend to those with one heavy-quark held (for bottom) and a Dirac held 
(for charm). As shown in Ref. [18], 

lim Zy (w) = Zy , 

mcCL^O " '' 

lim Zy^{w) = Zy^, 

rricCi^O 

\im Zv^{w)C]f^{w) = Cy {w)- 

mca^O ^ 
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(B33) 

(B34) 

(B35) 











the unbarred coefficients have no w dependence [27]. In practice, the error in these equations 
is of order as{a)mf.a. We shall neglect this contribution and use 


Pv^jw) _ Pv^, (w) _ 
Pv4w) ~ PV||H ■ 


(B36) 


The one-loop calculation of the right-hand side can be done at zero recoil and is, thus, much 
simpler. The one-loop result is given in the right-most column of Table V. To account for 
the error due to the neglected O^agrrica) contribution, as in Eq. (B30) we consider the size 
of the one-loop coefficient for the range of 6-quark masses used in this calculation. Ending 
< 0.352. With pmLx = 0.352 we take the error as 


± ay(2/a)p|iL”^2ca • 


(B37) 


3. Discretization errors at nonzero recoil {w > 1) 


Power-law discretization efiects arise from the higher-dimension terms in Eqs. (B20) 
and (B21). The discretization errors can be found by comparing the HQET description of 
lattice gauge theory to that of continuum QCD, as follows: substitute Eqs. (B18) and (B19) 
into Eqs. (B16) and (B17), multiply by the matching factors as in Eqs. (2.19) and (2.20), 
and form the combinations in Eqs. (2.22) and (2.23). The resulting HQET descriptions of 
the form factors are 


_ on _ 1 _ _ _ _ 

h+{w) = C+{w)E{w) + (nr (ic) [26(w;) + Ae(«^)] + , (B38) 

h_{w) = l{w + l)C_{w)E{w) + [2e3(w) - M{w)] - 

(B39) 


where^ 

By [w) 


B 


(/)mis, 


W) 


+ {w + l)CLGT(^) 2 j ’ 

( C^y^{w)Cy^,{w)\ 


(B40) 


B 


(/)LGT 


W) 


+ {w + l)C^A^{w) 

1 


C+{w) — wC_{w) -I- 


w-lC^^'^{w)CvAw) 


w 


dIOlgT/ n 

BV (w) 


4 (tc -f l)C'y'^'’"(M;) 


CJw) - 


2 CLGT(^) 

CLGT(^) 


(B41) 


As long as the matching of the dimension-three currents is carried out to order af, the parts 
of Eqs. (B40) and (B41) entailing the C coefficients collapses such that 


^Wmis ^ ^ (9('c^min(fc,r)+l^^ 


(B42) 


^ The continuum QCD analogs of Eqs. (B40) and (B41) can be obtained by erasing the superscript “LGT” 
and simplifying with Eqs. (B20) and (B21). The result becomes, as expected, trivial. 
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where the dimension-four currents have been matched through order In particular at 
the tree level {k = 0), 


^mis 

mis 


1 1 

-it-, 

2mc3 2mb3 

0 , 


(B43) 

(B44) 


while in continuum QCD, B± = l/2mc ± l/2mf, and B'^ = 0. Thus, we have tree-level 
matching in the dimension-four currents, with errors from this source of the form 


error 3 ,± = [/ 3 (moca) ± h{mQba)] Aa. (B45) 

Here 0 / 3 ( 77100 ) = l/ 2 m 3 — 1/2777.2, and the factor of A is a power-counting estimate of the 
HQET matrix element; A is the scale of nonperturbative QCD as it pertains to heavy-light 
mesons, roughly the difference between the heavy-light-meson and heavy-quark masses. 

Another discretization error arises from the function ^(tu) in lattice gauge theory and 
continuum QCD. In LGT, the kinetic and chromomagnetic masses appear. In this way, one 
hnds that the mismatch in in H(w) yields an error 

errors = [fsirrioca) /^(TTrofeo)] {w - l)Ao, (B46) 

taking the functions Ai to be of order A and building in the fact that the contribution 
vanishes as w —)■ 1. Similarly to above, a/s( 777 oa) = l/2mB — 1 / 27772 , which, for our choice 
of cswj is of order 0 ;^. 

Combining the two kinds of errors (© means to add in quadrature), 

h+{2.22) — / 7 +(cont.) = errors © — l)error 3 ^_|_, (B47) 

h_(2.23) — h_(cont.) = error 3 ^_. (B48) 


Because C- vanishes at the tree level, the contribution to the error in h- from C-errors is 
suppressed by an addition factor of and, thus, omitted here. Note that error 3 ^+ in h+{w) 
is multiplied by (tn — 1), whereas error 3 _ in h-{w) is not; cf. Eqs. (B38) and (B39). Our 
choices for the functions /s(^oa) and f^imoa) are discussed below; cf. Eqs. (B58) and (B59). 


4. Discretization errors at zero recoil (rc = 1) 


Because the next-to-leading-dimension discretization effects are suppressed by the 
next-to-next-to-leading-dimension effects may be of the same size. This is especially true 
at zero recoil, where the next-to-leading contributions to vanish. To capture the leading 
discretization errors of h_|_(l), therefore, one needs the dimension-five temporal vector current 
(with v' = v) [39]: 


^^— 3 

ry Tr4 y T / — © ^ h \ ( 1 , 1 ) ^-L^V (1,1) Cy D J_S ajS D jby 

Zy,V - -Zy^v ■ r - 2m3, 2m,, + 2™* 


(B49) 


( 0 , 2 ) ( 0 , 2 ) ' Bby ( 0 , 2 ) Cyi^by 

' ^ 87772 ^, + 


aEb 


^D\b 

( 2 , 0 ) CyD^b^ ( 2 , 0 ) CyS ■ Bby ( 2 , 0 ) Cyilpb 

' sA^ + AAA 


olEc 
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and similarly for —v ■ V. For the currents defined in Sec. II, as well as for the continuum 
currents, the r^-coefficients and ^-coefficients in Eq. (B49) all take the form 1 -|- 0{as)- The 
r^-like coefficients and associated masses in Eq. (B49) drop out of the analysis. 

From Eqs. (7.23)-(7.29) of Ref. [39], the HQET expansions through 0{]s?) of the matrix 
elements are 

= + (B50) 

where rjv is an HQET-to-QCD matching factor that starts with 1 in perturbative QCD, and 


'kS’ = 1 - - BA-AbE - + 3R2), 




- -I 

“ 2^3 


(1,1) 2 (1,1) 2 ' 

Vtt - 4s Fg 


(B51) 

(B52) 


where D, E, Ri, R 2 , and jjq are HQET matrix elements of order Also, 


A, 


1 

2 mic 


1 

2 mib' 


I 


2,5,3, 


(B53) 


are combinations of the mass coefficients in Eqs. (Bl) and (B14). Beyond the leading 1, the 
terms in VFoo^ come from double insertions of the kinetic and chromomagnetic interactions. 
Wqq stems from the dimension-Eve currents in Eq. (B49). 

Taking the difference between these expressions and the analogous ones for continuum 
QCD, one sees that the error in kFog^ stems from 


1 Zb 

2mBh 2m2h 


afBirnohO). 


(B54) 


The coefficients = 1 -|- 0{as)] also l/m^h + 0{asa) [compare Eqs. (B12) 

-( 2 ) 

and (B13)]. Thus, the error entering IFoo stems from 

A^yV^ - A4{/4^ = 2a[/3(moca) - /3(moba)]A2, (B55) 

with /3 of order for our choices. Thus, errors in pv^y/R^ stem from the mismatches 


IEo?^(LGT) - IFoo^(cont.) = -aA 2 [4(mo,a) - /^(moba)] {R, + 3^2 + 35), (B56) 

IfS^LGT) - lTi?(cont.) = -aA 2 [fsimoca) - fsimoba)] {pi - pi). (B57) 


In estimating heavy-quark discretization errors, we use these results at tc = 1, where the 
more generic effects in Eqs. (B47) and (B48) are much smaller. 


5. Numerical estimates 

For the mismatch functions Jb and in Eqs. (B46) and (B45), and in Eqs. (B56) 
and (B57), we use the functional forms [23] 



(B58) 

fRrriQa) = ^- - -r. 

^ ^ 2(2 + moa) 

(B59) 


® Ref. [39] used a notation setting = —Ai and fie = 3A2. 
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TABLE XIII. Absolute difference of h±{w) from mismatches in the heavy-quark Lagrangian and 
current. We take A = 450 MeV, = 0.424 GeV^, and = 0.364 GeV^. We further estimate 
the quantity |i2i -|- 3i?2 + 3E| with A^. The totals are obtained from Eqs. (B47), (B48), and 
(B56) and (B57) for h+{w), h-, and /i+(l), respectively. The column for /i+(tc) must be multiplied 
by {w — 1). The difference is estimated using the a = 0.09 fm lattice as a baseline. 


a (fm) 

av{q*) 

UlQb^ 


h+iw) 

1 

< 

/i+(l) 

0.120 

0.300 

2.462 

0.532 

-0.0095 

-0.0030 

-0.0011 

0.090 

0.261 

1.664 

0.362 

0.0000 

0.0000 

0.0000 

0.060 

0.220 

1.123 

0.240 

0.0109 

0.0021 

0.0011 

0.045 

0.198 

0.808 

0.176 

0.0160 

0.0029 

0.0016 


TABLE XIV. Absolute error on h±{w) from mismatches in the heavy-qnark Lagrangian and cur¬ 
rent. We take A = 450 MeV, = 0.424 GeV^, and /Jq = 0.364 GeV^. We further estimate 
the quantity |i?i -|- 3i?2 + 3E| with A^. The colnmns for h+{w) correspond to the chromomagnetic 
mismatch [“B”, Eq. (B46)], the current mismatch [“3”, Eq. (B45)], and their quadrature sum 
Eq. (B47)]; these columns must be multiplied by {w — 1). The column for h-{w) comes from the 
mismatch in Eq. (B48). The columns for /i_|_(l) correspond to the second-order mismatch of the 
Lagrangian [“VEqq^”, Eq. (B56)] and the second-order mismatch of the current [“VEog^”, Eq. (B57)], 
and their quadrature sum (“©”)• 


a (fm) av{q*) 

TTT-ofoft 

mocO 

h+{w) 

B 3 © 

h-{w) 

\/w 

h+{l) 
wg © 

0.120 

0.300 

2.462 

0.532 

0.0382 0.0125 0.0402 

0.0069 

0.0033 0.0005 0.0033 

0.090 

0.261 

1.664 

0.362 

0.0293 0.0092 0.0307 

0.0040 

0.0023 0.0003 0.0023 

0.060 

0.220 

1.123 

0.240 

0.0190 0.0057 0.0198 

0.0019 

0.0012 0.0001 0.0012 

0.045 

0.198 

0.808 

0.176 

0.0141 0.0041 0.0147 

0.0010 

0.0007 0.0001 0.0007 


To estimate the HQET matrix elements, we take A = 450 MeV,® 

= |(M^* - M|) = 0.364 GeV^ = (603 MeV)^ (BOO) 

/i^ (1 GeV) = 0.424 ± 0.042 GeV^ = (651 ± 32 MeV)A (B61) 

We do not have estimates for D, E, Ri, and R 2 as good as Eqs. (B60) and (B61), but in 
Ref. [5] we found that we could explain the discretization effects at zero recoil in B ^ D* 
with |i?i -b 3 R 2 + SE\ < (450 MeV)^. 

We take the typical av{q*) to be 0.262 on the a ~ 0.09 fm lattices, and we use one-loop 
running to obtain av{q*) at the other lattice spacings. 

In Table XIII, we show results from using these inputs to compute the differences that 
Eqs. (B56) and (B57) predict, using the lattice with a ~ 0.09 fm as the baseline. The 
estimates of the differences are compatible with the lattice-spacing dependence that can 
been seen for tc = 1 in Fig. 7, and can be inferred for w > 1 from Fig. 3. For example, 
the error in h+{w) grows slowly with w, both from Table XIII (adding in quadrature the 


® Here, 450 MeV is not an estimate of Mb —mb, but simply a practical number for power-counting estimates. 
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right-most column with (tc — 1) times the hfth column) and Fig. 3. Because the differences 
from lattice to lattice are well described by the theory, we can proceed to use the same 
ideas to estimate the difference from each lattice to the continuum. The results of these 
calculations are shown in Table XIV. For our hnal error estimates on the vector and scalar 
form factors, we take the absolute errors on and in Table XIV at a ~ 0.06 fm, and 
combine them in quadrature following Eqs. (2.5) and (2.6) that relate /+ and /o to and 
h_. The resulting expressions for the absolute errors as a function of recoil are 


error I = 


erroro = 


0.0198(u;- 1) 


0.0198(u;2 - 1) 


1 -I- r 
^/f 


1 + r 


© 

1 — r 
0.0019 ^ 

27^ J 

- © 
r 

0.0019(m; - 


© 0.0012 


1 + r 


a/t 


© 0.0012(u; + 1) 


a/t 
1 + r ’ 


(B62) 

(B63) 


These lead to estimates that range from 0.1-0.4% for both /+ and /o in our range of simulated 
lattice w values. 
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